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' Abstract 
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^ ! The dipole subtraction method for calculating next-to-leading order corrections 

I in QCD was originally only formulated for massless partons. In this paper 

we extend its definition to include massive partons, namely quarks, squarks 
and gluinos. We pay particular attention to the quasi-coUinear region, which 
gives rise to terms that are enhanced by logarithms of the parton masses, M. 
By ensuring that our subtraction cross section matches the exact real cross 
section in all quasi-coUinear regions we achieve uniform convergence both for 
hard scales Q ~ M and Q 3> M. Moreover, taking the masses to zero, we 
exactly reproduce the previously-calculated massless results. We give all the 
analytical formulae necessary to construct a numerical program to evaluate the 
next-to-leading order QCD corrections to arbitrary observables in an arbitrary 
process. 
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1 Introduction 



Hard-scattering production of heavy quarks and strongly-interacting heavy particles (such 
as squarks, gluinos and so forth) is of topical interest for Standard Model (SM) and beyond 
Standard Model physics at high-energy colliders (see, e.g. Refs. [|l], ^ and references therein). 
In order to have reliable theoretical predictions, it is important to control and explicitly 
compute QCD radiative corrections to these production processes beyond the leading order 
(LO) approximation in perturbation theory. These computations are very involved, as 
has been known since the first next-to-leading order (NLO) calculations of heavy-quark 
production in hadron collisions [0, ||. 

In the case of hard-scattering processes that involve only massless QCD partons, the 
practical feasibility of higher-order computations has been highly simplified by the develop- 
ment of general algorithms ^, 0, ^, |^, |TU|, |TT|, |T^] to perform NLO calculations. These al- 
gorithms, based either on the phase-space slicing method |jl3| or on the subtraction method 
|T^, start from the process-dependent QCD matrix elements and apply process-independent 
procedures to isolate and cancel the infrared (soft and coUinear) divergences that appear at 
intermediate steps of the calculation. The final output is a process-independent recipe to 
construct a modified version of the original matrix elements. The modified matrix elements 
can directly be integrated (usually by numerical Monte Carlo techniques) over the relevant 
and process-dependent phase space to compute any infrared and collinear safe observable 
at NLO. 



The NLO algorithm based on the dipole formalism was fully worked out in Ref. . 
In recent years, it has been implemented in general purpose Monte Carlo codes for the 
calculation of 3-jet and 4-jet observables |jT5|, in e+e~ annihilation, 2-jet [^, |T7| and 



3-jet JT^ final states in deep inelastic lepton-hadron scattering, production of three jets [p!8[] , 
vector-boson pairs [1^, vector boson plus massless hh |20| and colourless supersymmetric 
particles |21[] at hadron colliders. It has also been applied to the computation of QCD 
corrections to specific processes, such as 4-fermion final states and forward-backward 
asymmetries [EH] at high-energy e^e~ colliders and to radiative quarkonium decays 



Recent activity has been devoted to extending the dipole formalism by including the 
effect of massive partons. A first step was performed in Ref. p5|] (see also Ref. p6| for the 
particular case of small fermion masses) by developing the dipole formalism to compute 
NLO QED radiative corrections to electroweak processes and thus by considering photon 
radiation from massive charged fermions in arbitrary helicity configurations. The QED 
version of the dipole formalism was applied to study single-photon radiation in the scatter- 
ing processes 77 ff, e~7 e^7 and ^^fi^ — > z/gZ/e [11], to evaluate NLO electroweak 
corrections to e+e~ WW 4/ |^ and to W^-boson production at hadron colliders 
28| and to compute NLO QCD corrections to e~^e~ — >■ Htt (and Hbb) in the SM and its 



minimal supersymmetric extension |3^. A QCD extension of the results of Ref. |rT| to pro 



cesses involving heavy fermions has been presented in Ref. In this paper we present 
our version of the extension of the massless dipole formalism to QCD heavy partons. A 
small part of the results were anticipated in Ref. Full results of our formalism have 
already been applied to the computation of the NLO QCD corrections to the associated 



production of the SM Higgs boson and a tt pair in hadron collisions [32 
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In extending the NLO dipole formalism |11| to arbitrary processes with massive par 



tons, we devote particular attention to its behaviour in the massless limit. As long as 
we are interested in isolating and cancelling the infrared divergences in the cross section, 
the extension from massless partons to heavy partons mainly involve kinematical compli- 
cations due to the finite value of the parton masses. QCD radiation from heavy partons, 
however, can lead to contributions that, though infrared finite, are proportional to powers 
of InQ^/M^, where M is the parton mass and Q is the typical scale of the hard-scattering 
process. In kinematical configurations where Q ^ M, these logarithmically- enhanced con- 



tributions (as well as other types of constant contributions discussed in Sect. [21] ) become 
large and can spoil the numerical convergence of the calculation. In our formulation of the 
dipole method, we are not mainly concerned with the evaluation of these terms at NLO (or 
with their resummation to all orders), but rather in minimizing the instabilities that these 
terms can produce. For instance, in the case of cross sections that are infrared and collinear 
safe in the massless limit, these logarithmic terms cancel in the final NLO result, but they 
can still appear at intermediate steps of the calculation (e.g. in the separate evaluation of 
the real and virtual contributions) thus leading to reduced convergence in numerical im- 
plementations. These contributions can be singled out in a process-independent manner, 
since they are related to the singular behaviour of the QCD matrix elements in the limit 
M — >■ 0. This singular behaviour is controlled by quasi- collinear factorization formulae, in 
the same way that the infrared divergences are controlled by soft and collinear factoriza- 
tion formulae (see Sect. We exploit these factorization properties not only to cancel 
the infrared divergences, but also to stabilize the NLO algorithm with respect to large 
numerical contributions that can arise when the hard-scattering kinematics vary from the 
regions where Q ~ M to those where Q ^ M. In the case of processes with no initial-state 
hadrons (see Sect. [2?T| ), we set up the method such that, for those cross sections that are 
infrared and collinear safe in the limit of vanishing parton masses, the NLO algorithm is 
smooth in the massless limit and, more importantly, numerically stable for any values of 
the hard-scattering scale. In the case of hadron collision processes (see Sect. p.2|) , similar 
features are achieved provided the NLO partonic calculation is properly matched with the 
definition of the parton distributions of the colliding hadrons. 

In summary, we present a general extension of the dipole formalism that includes differ- 
ent species of strongly interacting massive fermions (quarks, gluinos) and scalars (squarks) 
with equal or unequal masses (see Sect. |^ and Appendix y). We develop the formalism by 
explicitly carrying out the analytical part of the NLO calculation for arbitrary infrared and 
collinear safe observables in lepton and hadron collisions (see Sect. The discussion of 
cross sections that involve fragmentation functions of massless or light (with respect to the 
hard-scattering scale) partons is only sketched, since it does not involve major technical 



complications with respect to the results presented here and in Ref. |TT|]. An important 
general feature of our treatment regards the identification and control of logarithmically- 
enhanced contributions of the type InQ^/M^, which can produce numerical instabilities in 
kinematical regimes where the hard-scattering scale Q is much larger than the value M of 
the mass of one or more heavy partons. 

This feature, which was also implemented in the case of photon radiation from mas- 



sive fermions |25], is the main overall difference with respect to other process-independent 



treatments p5| , ^ of massive partons in NLO calculations. There are, however. 



other general differences. For instance, in Refs. p5|, Q the soft divergences are regular- 
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ized by introducing an infinitesimal photon mass (a widespread practice in calculations of 
electroweak radiative corrections), while we use dimensional regularization. The authors 



of Ref. 1^ consider only fermions as heavy partons and do not consider the cases with 
different species of massive fermions of unequal masses. The general treatment of Ref. [0 
is based on the phase-space slicing method. The authors of Ref. pointed out the rele- 
vance of the quasi-collinear limit to control the large terms InQ^/M^ in the region where 
Q ^ M. They also discussed how these terms can be incorporated in the definition of 
parton distributions and fragmentation functions of heavy quarks but postponed a detailed 
treatment to future work. 

In the case of massless partons, the dipole formalism was fully developed and described 



in Ref. In the present extension to include massive partons, we closely follow the 



presentation in Ref. |[TT|]. Therefore, we do not repeat all the steps with the same amount 
of details as in Ref. |[Tl|. The outline of the paper is as follows. In Sect. ^ we discuss 
the general features of our implementation of massive partons in the dipole formalism. 
In Sect. ^ we recall the notation and its extension to the massive case. In Sect. ^ we 
discuss the singular behaviour of the real-emission matrix elements in the soft, collinear and 
quasi-collinear limits. In Sect. |^ we give the definition of the dipole factorization formulae 
that smoothly interpolate these limiting regions for any values of the parton masses. The 
proof that the definition consistently matches the singular behaviour of the QCD matrix 



elements is not explicitly given: it proceeds along the same lines as in Ref. [|TT|. In the case 



of vanishing parton masses, all the dipole factors reduce to those defined in Ref. |]TT|. In 
Sect, ^we also derive the factorization properties of the dipole phase space and perform the 
corresponding integration of the dipole factorization formulae. Sect. ^ contains the explicit 
expressions necessary for calculating QCD cross sections in different classes of processes at 
NLO accuracy. In Sect. |6.1| , we write down the final formulae needed to implement our 
method, while in Sects. |6.2| - |6.4| we sketch their derivation for processes with no (Sect. |6.2|) , 
one (Sect. |0| ) and two (Sect. |6.4|) initial-state hadrons. After reading Sect. |, a reader who 



is familiar with the method in the massless case and is interested only in setting up a NLO 



Monte Carlo program for a specific process can find all the relevant results in Sect. |6.1| and 
in one of the following subsections (depending on the specific process). Sect. |^ contains our 
conclusions. 

We leave some technical details and discussion of explicit examples to the appendices. 
In Appendix ^ we spell out the explicit expression for the integral of the eikonal term. 
Appendix ^ discusses a subtle point in the definition of the x-distributions related to the 
treatment of initial-state singularities in NLO calculations. In Appendix 0, we give the 
dipole splitting functions that are relevant for SUSY QCD calculations. Finally, we present 
the three simplest examples of our method in Appendix 0. 
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2 The general method and its features in the massless 
limit 



In this section we briefly review]^ the method of Refs. [|T^, |Tl|] and describe the main features 
of our extension to processes with massive partons. 



2.1 Cross sections without identified or initial-state hadrons and 
their smooth massless limit 

We first consider processes with no initial-state colliding hadrons and no tagged hadrons 
in the final state. This is the case, for example, of jet production in e~^e~ collisions. Our 
aim is to calculate a generic infrared- and collinear-safe cross section a with NLO accuracy. 
Assuming that at the LO there are m QCD partons in the final state, we schematically 
write 

' da^ + a^L", {2.1] 



^ = ^LO , ^NLO ^ / ^^B , „NLO 



where the LO contribution da^ is the fully differential Born cross section and the NLO 
correction a^^'^ comprises the real and virtual (one- loop) contributions dcx^ and dcx^: 



a^^^ = / da^+ da^ . (2.2) 

Jm+l Jm 

Here the notation for the integrals indicates that the real contribution involves m + 1 final- 
state partons (one QCD parton more than in LO), while the virtual contribution has the 
m-parton kinematics. 

After renormalization of the one-loop matrix element involved in da^, all the contri- 
butions to a are ultraviolet (UV) finite. The Born contribution dcx^ is integrable over the 
infrared (IR) region of the phase space, but the real and virtual contributions to a^^*^ are 
separately affected by IR divergences produced by soft and collinear partons. Although the 
IR divergences cancel in the sum on the right-hand side of Eq. ( p.2|) , the separate pieces 
have to be regularized before any numerical calculation can be attempted. Analytic con- 
tinuation of the integrals to c? = 4 — 2e space-time dimensions is the only known gauge- 
and Lorentz-invariant regularization procedure, but it prevents a straightforward numerical 
integration. Using dimensional regularization, the IR divergences are replaced by double 
and single poles, and 1/e, that have to be extracted and cancelled analytically before 
performing the limit e 0. 



The formalism of Ref. ITI] deals with the e poles by using the subtraction method. The 
general idea of the subtraction method is to introduce an auxiliary cross section do"^ that 
has the same pointwise singular behaviour (in d dimensions!) as dcr^. Moreover, do"^ has 
to be chosen simple enough, such that it is analytically integrable in d dimensions over 
the one-parton subspaces that cause the soft and collinear divergences. Thus, without 
performing any approximations, da^ is subtracted from the real contribution and added 
back to the virtual contribution. Since dcr^ acts as a local counterterm for dcr^, the 



§A concise overview of the dipole formalism can also be found in Rcf. 
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difference [da^ — da^] is integrable over the entire (m + l)-particle phase space in any 
number of dimensions and we can safely take the hmit e — > 0. Moreover, since the analytical 
expression of the performed integral over the singular subspace, da^, explicitly contains 
all the e poles that cancel those of the virtual term da^ , the sum [dcr^ + dcr^] is also 
(numerically) integrable in c? = 4 dimensions over the remaining m-parton phase space. 
The final structure of the NLO calculation is 



NLO 



m+1 



+ 



da^+ l^da' 



e=0 



(2.3) 



and the two terms on the right-hand side are separately integrable. Usually the integrations 
are not feasible in analytic form, but they can always be carried out numerically. For 
instance, the calculation can be implemented in a 'partonic Monte Carlo' program, which 
generates appropriately weighted partonic events with m + 1 final-state partons and events 
with m final-state partons. 

The method described so far is equally applicable to QCD calculations that involve 
massless and massive QCD partons. In principle, the extension from the massless to the 
massive case can be performed in a straightforward way. It is sufficient to overcome the 
(non-trivial) technical difficulties related to the generalization of da^ to massive partons 
and, in particular, to the analytic evaluation of the integral Ji da^ over the one-particle 
phase space with mass constraints. 



However, in our extension of the formalism of Ref. IT^], we require additional properties. 
To explain the reasons for these additional requirements, it is important to recall that the 
finite mass M of the QCD partons plays a different physical role in different physical 
processes. In some processes (e.g. the total cross section for heavy-quark hadropro duct ion 
0) the finite (and large) value of M has the essential role of setting the hard scale of the 
cross section. In these cases the massless limit M — > is IR unstable and the corresponding 
cross section cannot be computed in QCD perturbation theory. In other processes (e.g. the 
production of heavy-flavoured jets in e~^e~ annihilation [0]), instead, the hard scale Q 
is independent of the mass M and the latter has only the role of an auxiliary (though 
important) kinematical scale. These processes are IR stable in the massless limit, that is, 
when M — >■ the cross section is still infrared- and collinear-safe and, thus, perturbatively 
computable. 

The processes that are perturbatively stable in the massless limit are often studied in 
kinematical regions where the typical hard scale Q is much larger than the mass M of one 
(or more) of the heavy partons. In this regime the integral of the real term da^{M) of the 
NLO cross section in Eq. (|2.2|) leads to contributions of the type 



dcr^(M) 



m+1 



'dqi (qi)" 



1 



In -51- + 0(e) 



and 



m+1 



da^(M) 



(2.4) 



(2.5) 



where q_L generically denotes the typical transverse momentum of the heavy parton with 
mass M. Since these contributions are finite when e — 0, naively, they would not require 
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any special treatments within the subtraction method. However, this could lead to serious 
numerical problems. 



The problems are evident in the case of the contribution in Eq. ( p. 41) , which is very 
large when Q ^ M. When computing the NLO cross section, this large contribution would 
appear in the first term (the (m + l)-parton integral) on the right-hand side of Eq. (p.3| ) 
and it would be compensated by an equally large (but with opposite sign) logarithmic 
contribution arising from the second term (the m-parton integral). Owing to the presence 
of several large (although compensating) contributions, a similar naive procedure would lead 
to instabilities in any numerical implementations of the NLO calculation. The numerical 
instabilities would increase by increasing the ratio Q/M and, in particular, they would 
prevent from performing the massless limit. 



The contribution in Eq. (2^) may appear harmless, since it approaches a constant 
(finite) value when M/Q —>■ 0. However, the constant behaviour is obtained by combining 
the factor from the numerator of the integrand with the factor from the integral 

of the denominator. Owing to the presence of a linearly divergent (in the limit M^/ Q"^ —>■ 0) 
integral, the contribution in Eq. ( p.5| ) cannot be evaluated numerically by using standard 
Monte Carlo techniques, since its variance increases linearly with Q^/M^. Finite integrals 
with infinite variance (e.g. integrands with square root singularities) often occur in NLO 
calculations: they can be treated within Monte Carlo methods by applying importance 
sampling procedures. In this respect, the term in Eq. ( p.5|) poses no additional conceptual 
problems. However, our main point is that the variance of the integral in Eq. ( |2.5| ) is 
not uniform in M: it varies from a finite value when M = 0{Q) to a divergent value 
when M/Q 0. The presence of these type of contributions may thus prevent from 
straightforwardly constructing a 'partonic Monte Carlo' program that performs equally 
well in the two kinematical regimes M = 0{Q) and M/Q <^ 1. 

To avoid these numerical problems in the calculation of cross sections that are IR stable 
in the massless limit, we set up our massive-parton formalism by choosing the auxiliary 
cross section da^{M) in such a way that the following property is fulfilled: 

fcI«[K(^'))„.- K(^''))„ol = I«[K(^' = °))„„- ('J"'!^-' =»))„„] ■ 

(2.6) 

Note that, to avoid the problems related to the large logarithmic contributions in Eq. (|2.4|) , 
it is sufficient to impose that the integral of the subtracted cross section on the left-hand 



side of Eq. ( |2.6|) is finite when M — > 0. Equation (|2.6| ) is, instead, a stronger constraint. It 
implies that, in the evaluation of the subtracted cross section, the massless limit (or, more 
generally, the limit M/Q —>■ 0) commutes with the (m + l)-parton integral. This guarantees 
that [d(T^(M) — d(T^(M)] does not contain integrands of the type in Eq. (|2] 



Once Eq. (|2.6|) is satisfied, the (m -f- l)-parton and m-parton contributions to the NLO 
cross section in Eq. ( |2.3| ) separately have smooth behaviour when M — > 0. This behaviour 
helps to compute the cross section numerically in the kinematical regimes where Q ^ M. 
It can also be used to easily check that, in the massless limit, the massive-parton calculation 
correctly agrees with the corresponding result obtained in the exactly massless case. 

The IR divergent contributions to the real cross section dcr^ are process independent. 
By this we mean that, in the soft and coUinear limits, da^ is given by the corresponding 
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(process-dependent) Born-level cross section dcr^ times universal (process-independent) 
singular factors. Owing to these soft and coUinear factorization properties it 
is possible to give general prescriptions for constructing the auxiliary cross section da^ in 
a process-independent manner. 

Within the dipole formalism da^ is constructed by a sum over different contribu- 



tions, named dipoles. Each dipole contribution describes soft and coUinear radiation from 
a pair of ordered partons. The first parton is called emitter and the second spectator, since 
only the kinematics of the former leads to the IR singularities. The dipole configurations 
can be thought of as being obtained by an effective two-step process: using the Born-level 
cross section, an m-parton configuration is first produced and the emitter and spectator 
are singled out in all possible ways; then the emitter decays into two partons and the 
spectator, which contains information on the colour and spin correlations of the real cross 
section dcr^, is used to balance momentum conservation. The auxiliary cross section da^ 
can symbolically be written as 



da^= Yl da^®dKiipoie, (2.7) 

dipoles 



where the dipole factors dVdipoie describe the two-parton decays of the emitters. These 
factors are universal and can be obtained from the QCD factorization formulae (including 
the associated colour and spin correlations, as denoted by the symbol 0) in the soft and 
coUinear limits. 



The product structure in Eq. ( |2.7| ) is made possible by the factorization of QCD ampli- 
tudes on soft and coUinear poles and by a suitable factorized definition of the phase space. 
It permits a factorizable mapping from the (m + l)-parton phase space to an m-parton 
subspace, identified by the partonic variables in da^, times a single-parton phase space, 
identified by the partonic variables in dVdipoie- The single-parton phase space is process- 
independent: it describes the two-parton decay of the dipole and embodies the kinematical 
dependence on the degrees of freedom that lead to the IR singularities. This mapping 
makes dVdipoie fully integrable analytically and in a process-independent manner. We can 
symbolically write: 

/ da^ = ^ / da^ ® / dl^dipoie = / fda^ ® l] , (2.8) 

Jm+l jjp^i^^ Jm Jl Jm^ 

where the universal factor / is defined by 

1= J dKlipole , (2.9) 

dipoles ^ 



and contains all the e poles that are necessary to cancel the (equal and with opposite sign) 
poles in da^ . As a byproduct of this cancellation mechanism, Eq. (|2.8| ) can also be used to 
indirectly derive explicit information on the e singularities in the virtual contribution. Since 
dcr^ and da^ are respectively obtained from one-loop and tree-level QCD amplitudes, the 
factorization structure on the right-hand side of Eq. ( ^^ implies that the IR divergences of 
the one-loop amplitudes can be obtained from the corresponding tree amplitudes in terms 
of a universal factorization formula P, p]] . 
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The structure of the final NLO result is given as follows in terms of two contributions, 
^NLO{m+i} ^NLO{m}^ with (m + l)-parton and m-parton kinematics, respectively, which 
are separately finite and integrable in four space-time dimensions: 



cr 



NLO 



^NLO{m+l} _|_ ^NLO{m} 



(2.10) 



m+l 



da 



R 



e=0 



dV^, 



dipolc 



1 dipoles 



e=0 



da^ + da^ 



e=0 



Equation ( p.lO[ ) represents the dipole formalism implementation of the general subtraction 
formula (2.31). 



The exphcit expressions of the dipole factors dVdipoie and / for all the cross sections 
with massless QCD partons were obtained in Ref. 0]. In this paper we present their 



generalization to the case of massive QCD partons in the final state. As in Ref. [|rT|, the 
dipole factors dVdipoie are constructed by starting from the QCD factorization formulae in 
the soft and collinear limits. Moreover, to implement the smoothness condition in Eq. ( |2.6| ), 
the dipole factors have to match the behaviour of the QCD matrix elements in the so-called 



quasi- collinear limit [^]. The collinear limit describes the splitting process of one massless 
parton in two massless partons when the relative transverse momentum of the latter 
vanishes. The quasi-collinear limit (see Sect. |4.2| ) is a generalization of the collinear limit 
to the splitting processes of massive partons. It is obtained by letting and the parton 
masses vanish uniformly, i.e. we consider the limit q_L,M — >■ at fixed ratio q^/M. In 
the quasi-collinear limit, the integrands of the type in Eqs. (|2.4| ) and ( p.5| ) are formally 
regarded as being as singular as those that are proportional to 1/qi and lead to the 
collinear divergences. Therefore, constructing the dipole factors in Eq. ( |2.7|) in such a way 
that they have the same quasi-collinear limit as the real cross section da^, our (subtracted) 
{m-\- l)-parton contribution, (jNLO{m+i}^ -gq ( p.lO| ) fulfils the smoothness condition in 
Eq. (|2.6|) . In addition, the dipole factors of the present paper exactly coincide with those of 
Ref. [O in the massless limit. We do that for the practical purpose of facilitating a direct 



comparison with the results in Ref. |11 



Analogously to the massless-parton case, the cancellation of the e poles in the m-parton 
contribution a^^o^'^^ to Eq. ( glOD can be exploited to derive a universal factorization for- 
mula that controls the IR divergences of any one-loop amplitude with massless and massive 
partons. Having fulfilled the smoothness condition ( |2.6| ) and performing the analytical in- 
tegration over the one-parton phase space in Eq. ( |2.9| ) uniformly in the parton masses, our 
be used to obtain such a factorization formula in a form that explicitly ex- 
hibits not only the e poles, but also (i) the logarithmic terms (of the type In^ M and In M) 
that become singular when M — s> and (ii) the constant (finite) terms that originate in 
the calculation of the one-loop amplitude from the non-commutativity of the limits M ^ 
and e — > 0. This result was anticipated in Ref. pTI . 



2.2 Cross sections in hadron collisions 



The general procedure we have outlined so far applies to all processes with no initial- 
state hadrons (for instance, jets and heavy-quark production in e~^e~ annihilation). In 
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lepton-hadron and hadron-hadron collision processes, the hadronic cross section is com- 
puted by convoluting the corresponding partonic cross section (^(xphad) with the (process- 
independent) parton distributions f{x) of the incoming hadrons. Here, phad and x gener- 
ically denote the momentum of the colliding hadron and the momentum fraction carried 
by the colliding parton, respectively. The presence of initial-state partons, carrying a well 
defined momentum, makes the partonic subprocess collinear unsafe, so the partonic cross 
section cannot naively be computed in QCD perturbation theory. However, provided the in- 
coming parton is massless, the universal factorization theorem |^ of collinear singularities 



holds, and perturbative QCD can still be applied. The initial-state collinear divergences 
of the partonic cross section can be factorized and reabsorbed in the definition of the non- 
perturbative parton distributions. This procedure can consistently be carried out at any 
perturbative order. In particular, at NLO the perturbatively-computable partonic cross 
section is given by 



where 



a^LO^p) ^ f ^^R(p) ^ f ^^V(p) ^ f ^^C(p) . (2.11) 

J m+l Jm Jm 

I d(T^(p) = /^dx / da^{xp)T{x). (2.12) 

Jm Jo Jm 

Comparing Eq. ( p.2| ) with Eq. ( |2.11| ), we see that the latter contains the additional collinear 
counterterm da'^{p), which arises from the redefinition of the parton distributions. As 
symbolically written in Eq. ( p. 12 ), the collinear counterterm is given by the convolution 



of the Born cross section with a process- independent factor T{x), which is divergent when 
e ^ 0. 

We have emphasized the fact that the partons that initiate the hard-scattering subpro- 
cess have to be massless. The reason is that in the presence of two or more incoming massive 



partons the Bloch-Nordsieck mechanism |^ of cancellation of soft divergences is violated 



in QCD [^. The violation is due to non-cancelled soft divergences that are proportional 
to the mass M of the incoming partons and are process dependent. Although the violation 
occurs starting from next-to- next-to-leading order (NNLO), its process dependence spoils 
the general^] validity of the factorization theorem ||3^. Since the factorization theorem 
cannot formally be used to perturbatively define and compute partonic cross sections with 
incoming massive partons, we do not consider these cross sections. 

Note, however, that our distinction between massless and massive incoming partons 
has only a formal meaning. By 'mass' of the incoming parton we do not mean the actual 
physical mass of the parton. We only mean that the perturbative calculation is performed 
by setting this mass equal to zero. The practical calculational procedure is as follows. If the 
physical mass of a parton is much smaller than the typical value Q of the hard-scattering 
scale in the process, this parton has to be treated as a massless parton and included both 
in the initial and in the final state. If the case is not so, the parton has to be treated as 
massive and not considered as an incoming parton. We discuss this point further at the 
end of this subsection. 

The subtraction method and the dipole formalism can be used to compute the NLO 
partonic cross section in Eq. ( p.ll| ). Using the dipole formalism, the auxiliary cross section 

^By general we mean to all perturbative orders and for arbitrary non-vanishing values of the mass of 
the incoming partons. 
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do"^ still has the form in Eq. ( p.7|) and the sum over the dipoles is extended [11] to include 
initial-state dipoles that describe the coUinear decays of any initial-state massless parton 
(the emitter) into two massless partons. Furthermore, in considering the final-state emit- 
ters, the sum over spectators has to be extended to include initial-state massless partons, 
leading to dipole kinematics that differ from those with a final-state spectator. In the case 
with only massless partons in the final state, the initial-state dipoles were worked out in 

0. 



In this paper we extend the results of Ref. |11| to include the case with massive 



Ref. 

partons in the final state. 

Independently of the fact that the final-state partons are massless or massive, the initial- 
state dipoles can be integrated over the coUinear region. Their integral cancels the diver- 
gences of the coUinear counterterm da^{xp), and the final result for NLO cross section in 
Eq. ( |2.11|) can be written symbolically as 



cr 



NLO 



0" 



NLO {m+1} 



NLO{m} 



m+1 



+ 



da^(p) + da^(p) 




dx a^^^^"'\x;xp) 



dipolo 



(2.13) 



da^{xp) ®{P + K) (x) 



The contributions cr^'^o {m+i} ^^^^ ^NLO{m}^p^ (with (m + l)-parton and m-parton kine- 
matics, respectively) are analogous to those in Eq. (|2.10|) . The last term on the right-hand 
side of Eq. ( p.l3| ) is a finite (in four dimensions) remainder that is left after cancellation 
of the e poles of the coUinear counterterm in Eq. ( p.l2|) . This term involves a cross section 
^NLO {m} with m-parton kinematics and an additional one-dimensional integration 

with respect to the longitudinal momentum fraction x. This integration arises from the 
convolution of the Born- level cross section da^{xp) with x-dependent functions P and K 
that are universal and finite for e — >■ 0. The explicit expressions of P and K in the mass- 
less case were given in Ref. |]ll| . In this paper we obtain their generalization to the case of 



partonic cross sections with massive partons in the final state. 



In Ref. |]lT| an extensive discussion was devoted to a fully detailed treatment of cross 
sections with identified hadrons in the final state, i.e. cross sections that involve the intro- 
duction of parton fragmentation functions of the outgoing hadrons. For simplicity, in this 
paper we do not present an analogous discussion of these cross sections, since the exten- 
sion to include the associated production of final-state heavy partons is straightforward. 
More precisely, we can consider two main kinematical configurations according to whether 
the momentum Q of the final-state massive parton is a) not observed or comparable to 
its mass M and b) much larger than its mass. In the case a), the heavy parton does not 
require additional (with respect to other cross sections) special treatment, apart from it 
being included as spectator in the final-state dipoles that are introduced to describe the 
coUinear decay of the massless emitter that undergoes final-state fragmentation. In the case 
b), a proper treatment of the large contributions InQ^/M^ is required, which can be simply 
achieved by considering the massive parton as effectively massless and by introducing its 



(process-independent) perturbative fragmentation function |^ 



Considering the case of cross section calculations in processes with no initial-state 



hadrons, in Sect. 2.1 we have discussed how we set up the dipole formalism with mas- 
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sive partons to guarantee a smooth and numerically stable interpolation throughout the 
kinematical region where the hard-scattering scale Q varies from Q ^ M to Q ^ M. In 
hadron collisions, a QCD parton with a non- vanishing physical mass M has to be regarded 
as effectively massive or massless depending on whether Q ~ M or Q ^ M. Moreover, 
the transition between the former to the latter region needs the introduction and scheme- 
dependent definition of the parton distribution of the effectively massless parton. Evidently, 
a smooth interpolation between the two regions requires a proper and careful matching be- 
tween the definition of the partonic cross section and that of the parton distribution. This 
general theoretical issue has been the subject of many investigations in the recent litera- 
ture (see e.g. the list of references in Ref. [^), and several practical implementations are 
available, mostly in the context of deep-inelastic heavy-quark production. Once the match- 
ing conditions and the parton distributions within a certain scheme are precisely specified, 
smoothness conditions similar to Eq. ( ^.6| ) can be imposed also on NLO calculations of 
cross sections that involve parton distributions and fragmentation functions. The proce- 



dure outlined in Sect. ^TT] to guarantee the validity of Eq. ( p.6| ) can thus be applied to set 
up the dipole formalism such that it performs in a manner that is smooth and numerically 
stable with respect to large variations of the ratio Q/M. This is left to future studies. 



2.3 Summary of the general method 

The starting points of the calculation of QCD radiative corrections at NLO are the tree- 
level and one-loop matrix elements that respectively enter in the expression of the cross 
section contributions do"^ and da^ in Eq. ( |2.2| ) (or Eq. ( p.ll[ )). The separate integration of 
these matrix elements is not trivial because of their infrared divergences. The goal of the 
dipole formalism is to construct 'effective matrix elements' that can straightforwardly be 
integrated in four space-time dimensions. The final output of the formalism is summarized 
in Eq. (^TTUp (or Eq. ( |2.13| )), which give the expression of the NLO cross sections in terms of 
the 'effective matrix elements'. The formalism provides explicit expressions for the universal 
factors dVdipoie and I (and P, K), which we present in Sect. |[ Having these factors at 
our disposal, the only other ingredients necessary for the full NLO calculation are simply 
related to the evaluation of the original matrix elements. We need: 



• a set of independent colour projections of the matrix element squared at the Born 
level, summed over parton polarizations, in d dimensions; 

• the one-loop contribution da^ in d dimensions; 

• an additional projection of the Born-level matrix element over the helicity of each 
external gluon in four dimensions; 

• the real emission contribution dcr^ in four dimensions. 



Since the comment presented in the 'Note added' section of Ref. |Tl[] is valid independently 
of the value of the parton masses, this list can be simplified if one uses the dimensional- 
reduction scheme for regularizing the one-loop matrix elements. In this case the Born-level 
calculation in the first item above can directly be carried out in four space-time dimensions. 
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3 Notation 



3.1 Matrix elements 



We consider processes with coloured particles (partons) in the initial and final states. Any 
number of additional non-coloured particles is allowed, too, but they will be suppressed in 
the notation. Partons in the initial state are labelled by a, 6, . . ., partons in the final state 
by i,j, k, . . generically we write {i; a} for all partons. When we do not make a distinction 
between initial- and final-state partons, we use the common labels I, J, . . . and we shortly 
write {/} = {i; a}. 

The colour indices of the partons are denoted by c^, Cj, which range over 1, . . . , A^"^ — 1 
for gluons (or any other partons, such as gluinos, in the adjoint representation of the gauge 
group) and over 1, . . . , A^'c for quarks and antiquarks (or any other partons, such as squarks, 
in the fundamental representation). Spin indices are generically denoted by Sa, Si. As in 
Ref. [^, we formally introduce an orthogonal basis of unit vectors Ca}) ® \{si; Sa}) in 
the space of colour and spin, in such a way that an amplitude Ai^'^^''^^'^"''^''^ {{pf, Pa}) with 
definite colour, spin and momenta {pi,Pa} can be written as 



_^{c.s.;c.,s4(|p.pj) ^ [^n\A5^j ® {{s^■, Sa}\)\{t; a}) . (3.1) 

Thus \{r, a}) is an abstract vector in colour and spin space, and its normalization is fixed by 



including a factor of l/ync(6) for each initial-state parton b carrying ndb) colour degrees 
of freedom. Then the squared amplitude summed over colours and spins and averaged over 
initial-state colours is 

Y^^^\M{p^■,Pa})\' = {{^■,a}\{^■,a}). (3.2) 

Colour interactions at the QCD vertices are represented by associating colour charges 
Ti or Ta with the emission of a gluon from each parton i or a. The colour charge Tj = {T"} 
is a vector with respect to the colour indices n of the emitted gluon and an SU{Nc) matrix 
with respect to the colour indices of the parton i; analogously Ta describes gluon emission 
from the initial-state parton a. More precisely, for a final-state parton i the action onto 
the colour space is defined by 

(ci, . . . , Cj, . . . , Cm|T"|6i, . . . , 6i, . . . , 6m) = ^cifci • • • T^^bi ■ ■ ■ ^c^b^ , (3.3) 

where T^!^ is the colour-charge matrix in the representation of the final-state particle i, i.e. 
^cfe = i/cn6 if i is a gluon or a gluino, T"^ = t^^ if i is a (s) quark and T^^ = — if ^ is an 
anti(s)quark. The colour-charge operator of an initial-state parton a, for which a relation 
analogous to Eq. (|3.3| ) holds, is defined by crossing symmetry, that is by T^^ = —t^^ if a 
is a (s)quark and T^^^ = t^^ if a is an anti(s)quark. Using this notation, we also define the 
square of colour-correlated tree-amplitudes with m final-state partons, lAI^'^p, as follows 

\Mii\' ^ (jlMb^j m{{^,a}\T,■T,\{^,a})m 

= [MT-'--'"'-^-{{pf,Pa})]* T,^^a,Tb:a,MT-''^-^^-^-{{pf,Pa}) , (3.4) 
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and analogously for the cases in which j and/or k are replaced by initial-state partons. 

In our notation, each vector is a colour- singlet state, so colour conservation is 

simply 

(E^. + E^O |{^;4) = o, (3.5) 
^ j b ^ 

where the sum over j and h extends over all the external partons of the state vector a}) 
in the final and initial state, respectively. 

The colour-charge algebra for the product {Ti)^{TjY = Ti ■ Tj is: 

T,-Tj = T,-T, iity^j; Tl = C,, (3.6) 

and analogously for initial-state partons. Here Cj is the quadratic Casimir operator in the 
representation of particle i and we have Cp = Tj^^N^ — l)/Nc = {N^ — l)/{2Nc) in the 
fundamental and Ca = STrA^^c = Nc in the adjoint representation, i.e. we are using the 
customary normalization Tr = 1/2. 



3.2 Dimensional regularization, one-loop amplitudes and renor- 
malization 

We employ dimensional regularization in c? = 4 — 2e space-time dimensions to regulate both 
the IR and UV divergences. More precisely, we use conventional dimensional regularization 
(CDR), where quarks (spin-| Dirac fermions) possess 2 spin polarizations, gluons have d — 2 
helicity states and all particle momenta are taken as (i-dimensional. The li-dimensional 
phase space for m outgoing particles with momenta pi, . . . ,Pm, masses mi, . . . , rrim and 
total momentum P is denoted by 

i2nf5^''\p, + ...+Pm-P). (3.7) 



d4>m{pi, ...,Pm;P) 



n 



\ (2vr) 



Pi 

d-l 



The dimensional-regularization scale, which appears in the calculation of the matrix 
elements, is denoted by fi. The dependence on fi cancels after having combined the matrix 
elements in the NLO calculation of physical cross sections, although the latter eventually 
depend on the renormalization and factorization scales and np. To avoid a cumbersome 
notation, we therefore set fiR = /x, while /i and fip will differ in general. 

The virtual contribution do"^ to the NLO cross section is proportional to the real part 
of the interference between the tree-level and one-loop matrix elements. Although this 
interference is not positive definite, we denote it by |A^({pj,Pa})|(i_ioop)' since it represents 
the one-loop correction to the square of the tree-level matrix element. Actually, dcr^ is pro- 
portional to the renormalized one-loop correction |A^(fe,Pa})|(i_ioop)5 which is obtained 
from the corresponding bare quantity by adding ultraviolet counterterms that implement 
charge (coupling) and mass renormalization. The renormalized QCD coupling at the renor- 
malization scale fi is shortly denoted by ctg. It can be defined either in the customary MS 



renormalization scheme or in the renormalization scheme of Ref. 43 , which differs from 
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the MS scheme in the treatment of heavy partons and guarantees their decouphng in the 
infinite mass hmit. The renormahzed mass parameters (i.e. parton masses and related pa- 
rameters such as those which appear in Yukawa couphngs) are denoted by mi and can refer 
either to the MS renormahzation scheme or to the pole-mass definition. The results of this 
paper are independent of these different renormahzation prescriptions. In particular, since 
these prescriptions coincide in the limit of vanishing parton masses, they do not affect the 
massless limit of IR stable cross sections. 

Once renormahzation has been performed, the one-loop contribution |A^|^i_ioop) 
longer contains 1/e poles of UV origin, but it still contains 1/e poles and finite terms of IR 
origin. The latter depend on the dimensional regularization procedure used for evaluating 
the loop integral. The regularization-scheme dependence of the virtual corrections has to be 
consistently matched to that of the real ones, to guarantee the regularization-scheme inde- 
pendence of the NLO cross section. Since in this paper we deal with the real corrections by 
using the CDR scheme, we need the result for |A^|(i_ioop) within this regularization scheme. 
One-loop matrix elements are sometimes evaluated by using dimensional-regularization 
prescriptions that differ from CDR (the difference can be due to the dimensionality of the 
momenta of the external particles and/or to the number of polarizations of both external 
and internal particles). If |-^|^i_ioop) known in a scheme that is different from CDR, its 
expression in the CDR scheme can be obtained by introducing a correction proportional to 
the corresponding tree-level amplitude. The correction term for several different regular- 
ization schemes can be found in Refs. ?3| for the massless-parton case and in Ref. [ JT| 



for the general case with massless and massive partons. 



4 Factorization in the soft and (quasi- )collinear limits 

We consider a generic tree-level matrix element A4m+i with m + 1 massive or massless 
QCD partons in the final state and up to two massless initial-state partons. At NLO, 
the dependence of the squared matrix element |A^m+iP on the momenta of the final-state 
partons is singular in two different situations. The first situation occurs in the soft region, 
where the momentum of a final-state gluon tends to zero in any fixed direction. The second 
situation corresponds to the (quasi-) collinear region, where a final-state parton becomes 
collinear to another (initial- or final-state) parton. 



4.1 Soft limit 

In the soft region the momentum pj of a final-state gluon j tends to zero in any fixed 
direction q, i.e. pj = Xq with A ^ for fixed q. In this limit, the squared matrix element 
diverges as and its divergent part can be computed in terms of the eikonal current 



of the soft gluon As described in Ref. |[TT[] in more detail, using partial fractioning 



and colour conservation, the singular behaviour of can be written as a sum over 
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emitter-spectator (I-K) pairs, 



m+l,a... (• • • ; J; • • • ) Ct) • • • I I • • • ; J; • • • ) Ct) • • ■)m+l,a.. 
1 



(4.1^ 



A2 



m,a...{- ■ ■ 1 O'y 



\Tt-T 



K 



PiPk 



m, 



iPi+PK)q '2piq 



where | . . .)m,a... corresponds to the m-parton matrix element that resuhs from the original 
(m + l)-parton matrix element upon omitting the gluon j, and the indices /, K label both 
initial- and final-state partons. This formula is valid in the general case of massive or 
massless partons. 



4.2 Collinear and quasi- collinear limits 

In the massless case, the squared matrix element diverges when two external 

partons become collinear |3^. If at least one of the external partons is massive, the collinear 
divergence is screened by the finite value of the parton mass. Nonetheless, the cross section 
contribution of the matrix element from this phase-space region is strongly (logarithmically) 
enhanced when the parton mass becomes small. As discussed in Sect. |2]1|, to control these 
enhanced contributions we have to study the behaviour of the matrix element in the quasi- 
collinear limit |^. Since the quasi-collinear limit is a generalization of the customary 



collinear limit, the latter can in turn be obtained as a special case from the former. 

We consider two final-state partons i and j that can be produced by a QCD vertex 
through the splitting process ij i + j. The partons i and j have momenta pi and pj and 
masses rui and nij {pf = mf.p^j = m^). The mass of the parton ij is denoted by rriij and is 
constrained by rriij < "^i + f^ij- This constraint forbids the on-shell decay of the parton ij. 
We introduce the following Sudakov parametrization of the parton momenta: 



kl + z^ml -rnf nt" ^ _ ^1 + (1 - zfrnl - m) 



2pn ' ^ ' " ^ 1 — z 2pn 

(4.2) 

where p'^ is a time-like momentum (with = mfj) that points towards the collinear 
(forward) direction, n'^ is an auxiliary light-like vector (n^ = 0) and kj_ is the momentum 
component that is orthogonal to both p and n {pk± = nk± = 0). The invariant mass of the 
final-state partons is 

z{l — z) z 1 — z 

The quasi-collinear region is reached when k± becomes of 0(m) and small. This region 
can be identified by performing the uniform rescaling 

k± — s> Xk±, rrii Xrrii, mj Xrrij, rriij Xrriij (4-4) 

and studying the limit A — > 0. Neglecting terms that are less singular than 1/A^ the squared 



matrix element behaves as pi 



m+l,a... (• • • ; ^5 J; • • • ) • • • I I • • • ; ^5 j'j • • • ) C^) • • ■)m+l,a... 

1 Svr/i^^as 



A2 {pi+PjY -m^j 



2 m,a...{- . . , ij, . . . ; a, . . . i(2, {m}] e)| . . . . . . ; a, . . .)m,a...,(4.5) 
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where | . . .)m,a... corresponds to the m-parton matrix element that is obtained from the 
(m + l)-parton matrix element by replacing the parton pair z, j by the single parton ij. 

The kernel P~. ■{z,k^,{m]]e) in Eq. ( [4.5|) is the generalization of the d-dimensional 
Altarelli-Parisi splitting function from the coUinear to the quasi-collinear limits. As in 
the coUinear case, it depends on the momentum fraction z, on the transverse momentum fc^ 
and on the helicity of the parton ij in the m-parton matrix element. In the quasi-collinear 
case, it also depends on the masses mi,mj,mij of the partons involved in the splitting 
process. The mass dependence is shortly denoted by {m}. 



The generalized Altarelli-Parisi kernels for the QCD splitting processes Q 
Q ^ g + Q and g ^ Q + Q, are 



Q + g, 



{s\pQQ{z,k^,mQ;e)\s') 
{s\pQg{z, k±,mQ;e)\s') 
{fi\PgQ{z,k±,mQ;e)\iy) 



6ss' Cf 



1 + 
1 - z 

1 + (1 



— ez — 



g^^ - 4 



Aj _|_ 



+ 





(4.6) 


PQPg. 




' 


(4.7) 


P9PQ. 






(4.8) 



where s, s' and n, v denote the indices of the spin correlations for the parton zj, and the 
k^ and mass dependence of the scalar products 'p(Pj is given in Eq. (|4.3| ). Of course, the 
g ^ gg splitting function remains unaffected by the inclusion of mass terms. 



{^l\Pgg{z,k^■e)\v)=2C^_ 



-g 



1 



2(l-e)z(l-z) 



(4.9) 



We also need the spin-averaged form of the splitting functions. For quarks this is obtained 
by contraction with 5^8' /2 and for a gluon with on-shell momentum p, by contraction with 

1 , , . 1 



-d^uip) 



-gtiu + (gauge terms)] 



01-2^"'"^' 2(1 -e) 
where the gauge terms are proportional either to or to p'^ and 

- g^'d^^ip) =d-2, p^d^,{p) = 0. 
The spin-averaged QCD splitting functions only depend on z, e and the ratio 



mf + m| 



{Pi + PjY - w 



2 ' 



and are denoted by {P~.^{z; e; ^ifj)). They explicitly read 



{pQQ{z]e]^ilg)) 
{pQg{z;e;iilQ)) 
(^.Q(^;e)) 
{Pgai^;^)) 



Tr 
2Ca 



^-e{l-z)-2f,^ 
1 + (1- 



Qg 



z 

2z(l 



— ez — 2fi 



gQ 



z - 



z 



1 - z 
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+ 



+ z(l-z) 



(4.10) 
(4.11) 
(4.12) 

(4.13) 
(4.14) 

(4.15) 
(4.16) 



These generalized, spin-averaged splitting functions and some of their counterparts in SUSY 
QCD have already been given in Ref. |^T|. The four- dimensional QED analogue of these 



functions describes the splitting processes — -|- 7, 7 — >• e+ -|- e and is well known 



Setting rrii = rrij = rriij = in all the expressions of this subsection, we recover the 
known expressions for the collinear limit (see e.g. Ref. [pi]]). 



The matrix element Aim+i,a... diverges also when a massless final-state parton i becomes 
collinear to a massless initial-state parton a. The singularity is related to the splitting 
a ^ ai + i, where ai is the incoming parton of the related m-parton process. Since in this 
paper we restrict ourselves to massless incoming partons, none of the partons a, i, ai is 
massive and we can completely take over the factorization for massless partons described 



in Ref. There, the collinear limit is defined hj kj_ ^ in 



P'^=il-^)p'a+kl-^^, PaP^ = -,^^^y (4.17) 

1 — X IpaU 2(1 — X) 



In this limit |A^m+iP behaves as 



m+l,a... (...,i,...,fl,...||...,^,...,Cl,.. 



(...;a^,...|P ~.(x,fc^;e)|...;m,...) ~. , (4.18) 



where I . . .) ~ corresponds to the m-parton matrix element that is obtained from the 

I ' m,ai... ^ 

(m -|- l)-parton matrix element by replacing the parton pair a, z by the single incoming 
parton az, which has the reduced momentum xpa with < a; < 1. 



5 Dipole factorization formulae 



In the actual calculation of cross sections, Eqs. ( [4.1|) , ( [4.5|) and ([4.18|) cannot straightfor- 



wardly be used as 'true' factorization formulae, because the momenta of the partons in the 
matrix elements on the right-hand sides are unambiguously defined only in the strict soft 
and (quasi-) collinear limits. In this section we define dipole factorization formulae that, 
besides having the correct limiting behaviour, implement momentum conservation away 
from the limits. Momentum conservation is implemented not only exactly but also in a 
factorizable way. The dipole factorization formulae can thus be used to introduce the dipole 
factors dVdipoie that are needed to construct (see Eq. (|2.71) ) the auxiliary cross section dcr^ 
and its related integrated counterpart in Eq. 



We present the dipole factorization formulae as in Ref. |jTTI 



P = E E ^^.--^ + E E + E E + E E ^"'^ + • • • ^ (s-i) 

i,j k^ij i,j o, a,i j^i a,i bj^a 

where the terms on the right-hand side approximate the matrix element in different singular 
regions and the dots stand for terms that lead to finite integrals. Four different situations 
can be distinguished for the emitter-spectator pairs (see Fig. |lD, since either one can 
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Figure 1: Effective diagrams for the different emitter-spectator cases. 



belong to tlie final or initial states. The first two terms on the right-hand side control 
the singularities of the (m + l)-parton matrix element when two final-state partons i and 
j become (quasi-) collinear: the emitter is a final-state parton and the spectator can be 
either in the final {T>ij^k) or in the initial {l^lj) states. The third and fourth terms on the 
right-hand side control the singularities of the {m+ l)-parton matrix element when a final- 
state parton i and an initial-state parton a become collinear: the emitter is an initial-state 
parton and the spectator can be either in the final ("P^*) or in the initial [T>"'^'^) states. 
When the parton i is soft, all the four dipole functions V become singular. These four 
types of dipoles are considered in the following subsections. 

The dipole factors introduced in Ref. differ from ours in several respects. For 



instance, they do not deal with the case of massive partons with unequal masses and, 
typically, they treat the kinematic recoil differently. The main overall differences arise from 
the fact that they are not aimed to control the quasi- collinear region. 



5.1 Final-state emitter and final-state spectator 



The dipole contribution (see Fig. |l]) to the factorization formula ( |5.1|) is 

2~ "*(•••; iji . . . , k, . . . I ■ ■ ■ yij , ■ ■ ■ , k, . . .)fn ■ (5.2) 



{Pi+PjY -mj. T 
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The final-state parton momenta pj, pj andp^ have arbitrary masses and their total outgoing 
momentum is denoted by Q, 

pl = ml , p'^.=m], pl = ml, Q = Pi + Pj + Pk ■ (5.3) 
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The m-parton matrix elements in Eq. ( ^.21) are obtained from the original (m + l)-parton 
matrix element by replacing i and j with the parent parton ij in the splitting process 
ij ^ i + j and by changing the momentum of the spectator parton k. The on-shell mass 
of the emitter ij is rriij and we only consider situations with rriij < rrii + rrij . In Eq. ( |5.2| ) 
we have split off the colour structure from the spin functions ^ij,k, which are given in 



Sect. 5.1.2 



5.1.1 Kinematics and phase-space factorization 

The auxiliary momenta Pij and pk of the emitter and spectator are defined in such a way 
that Pij plays the role of p in the quasi-collinear limit ( [4.2[ ), but at the same time they obey 
their mass-shell conditions and total momentum conservation, 

pI = ml , Pl = ml, = pfj + p'^ . (5.4) 

For later use, we introduce the rescaled parton masses /x„ and the relative velocities Vp^g 
between two massive momenta p'^ and q'^, 



_ f^n { _ ■ ■ h, ■ ■\ _ 



1 - S^. (5.5) 



The relative velocities can also be written as 



in terms of the customary triangular function X{x, y, z), 

X{x,y,z)=x +y +z — 2xy — 2xz — 2yz . (5.7) 
For instance, the velocity Vij^k between pij and pk is given by 

= 1 -2 72- ■ (5-8) 

Our definition of the momenta Pij and Pk in terms of the original momenta Pi, pj and 
Pk (or Q =Pi+Pj +Pk) is: 

VA(Q^(p.+P.)^mDV'" J ^' 

P^=Q^'-P^ (5.9) 



and it coincides with the one given in Ref. [^] for the specific case with = rriij consid- 



ered there. The definition in Eq. (|5.9|) is symmetric with respect to i ^ j and is uniformly 
applicable to any configuration of the parton masses m„ (n = i,j,k,ij). Different defini- 
tions of the momenta pij and pk can be introduced. The authors of Ref. consider only 
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the specific cases with {nii = rrij 



rrii 



0,mk 7^ 0}, {iTLi = mij,mk = rrij = 0} and 



{rrii = rriij = ruk, rrij = 0}: in the last two cases their definition coincides with ours. 

Next we exactly factorize the three-particle phase space d(j){pi,pj,pk; Q) of Eq. ( p.7|) in 
terms of the two-particle phase space d(f){pij,pk; Q) and a single-particle phase-space factor 

[dpiipij.Pk)], 



d(l){pi,pj,pk; Q) = d(t){pij,pk] Q) [dpi{pij,Pk)\ 0(1 - /ij - /ij - /ifc) • 



(5.10) 



In App. B of Ref. p5[ the derivation of this phase-space splitting was outlined for the 
case with mi = rriij in four space-time dimensions. The generalization to rrii 7^ ruij in d 
dimensions is straightforward; we present only the results. We obtain 



:i - - - f^ir'- 



1 + 26 



|[dp,(p,,,p,)] =^(27r)-3+2^(g2)i- 

dyij,k (1 - [fJ'l + fJ'^j + (1 - Ati - Atj - fJ'l)yij,k 



X 



dzi [z+{yij^k) - Zi] ' [zi - z^{yij^k)y 



where the variables Zi and yij^k are defined as in Refs. [|ri|, [25|] : 

PiPk ^. PiPj 



Zi 1 Zj 



PiPk + PjPk 
Their integration boundary is given by 

2fiijij 



yij,k 



pipj + PiPk + pjpk 



(5.11) 



(5.12) 



y- 



y+ = i 



2 ' 



z±{yij,k 



2/i^ + (1 - /i- - /ij - ^4)yij,k 
2[/^i + + (1 - /i? - - 



-(1 ± Vij^iVij^k) , 



(5.13) 



where Vij^i ivij^k) is the relative velocity between pi +pj and Pi (pk). The relative velocities 
are functions of yij^k and explicitly read 



"^ijtk 



[2fil + (1 - /i^ - /ij _ ^2)(i _ ^^^.^)]2 _ 4^: 



(1 -/^f -/^i -^i)(l -l/y,fc) 

1 - ^2 _ 2 _ ^2)2^2 _ 4^2 2 



(5.14) 



The d'^"^!] integration extends over the solid angle perpendicular to pij and pk and thus 

271 



d'^-^Q 



vrTfl -e) 



(5.15) 
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5.1.2 The dipole splitting functions 



We give the functions Vj^ ^ in Eq. ( |5.2| ) for the three QCD splitting processes ij ^ i + j: 

• Q gipi) + Q{pj). mi = and mj = rriij = mq, 

• g —> Qipi) + Qipj): rrii = rrij = mq and rriij = 0, 

• 9 ^ g{Pi) + g{Pj)- = rrij = rriij = 0. 



The case Q — gQ is formally identical to Q — ^ gQ. The spectator mass m^, may be zero 
or non-zero in all cases. Note that in general Vj^ ^ is non-diagonal in the helicity space 
of the parton ij. For the analytical integration over the singular degrees of freedom (see 
Sect. |5.1.3|) we also need the spin-averaged functions (Vjj^fc)- Denoting s,s' or yU, z/ the 
helicities of ij = Q or ij = g in {. . . ,ij, . . .\ \ . . . ,ij, . . .), we define the dipole functions 
Vjj^fc as follows: 



{s\'VgQ,k\s') = M^i'^^a^C^ I 



(1 - yij,k) Vij^k 



PiPj 



(5.16) 



'^ij,k 



2k 



m 



1-e 



Z4-Z_ 



Q 



{Pi + Pi) 



{Pi + PjY 

2 



~(m) a ~(iri) a 

Zi Pi-z) 'p^j 



~(m) u ~('m) jy 
A Pi - Zj 'Pj 



,(5.17) 



1 f 2 

{^QQ,k) = STc^'^'a^Tji <^ 1 - -— 



Zi{l - Zi) - (1 - k)z+z. 



{lJ'\'ygg,kW) = ISvr/i 'OsCa { -g 



+ (1 - - fil)yij,k 
1 



(5.18) 



+ 



2 — KZ-i-Z- 



1 - Zi{l - yij^k) I - Zj{l - yij^k) Vij^k 



+ 



1 1-e 

Vij,k PiPj 



{^gg,k) = IQlTn'^^asCA | ^ _ ^ 



~(m) a ~{m) a 

A Pi-Zj 'Pj 
1 



■~(m) y ~(m) y 
Zi Pi - Zj 'Pj 



(5.19) 



+ 



yij^k) I - Zj{l - yij,k) 
Ziil - Zi) - (1 - k)z+z_ - 2 1 



Vij,k 



(5.20) 



where z± is given in Eq. ( |5.13| ) and the new variables z^"^^ and Zj^' are 

~(m) ~ 1/-. N ~(m) ~ 1/-. N 

Zi ' =Zi--{l- Vij^k) , Z) ' = Zj - -(1 - Vij^k) ■ 

The constant k is a free parameter, which only redistributes non-singular contributions 
between the different terms within square brackets in Eq. (p.3|). Choosing different values 



(5.21) 
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of K, we can simplify the expressions of either the subtraction term (k = 0) or its integral 
(k = 2/3). Of course, the final result for o"^^*"* must not depend on i.e. its k- independence 
can serve as a consistency check of the calculation. The velocity factors Vij^k/vij^k and 
1 / Vij^k are also harmless in the singular limits; they are introduced in the definition of the 
functions Vj^ jt only to simplify their integration in analytic form. The definition of the 
four- dimensional V^g ^ in Ref. differs from ours by some harmless velocity factors. 



In App. ly, we give the dipole splitting functions 'Vij,k that are relevant for SUSY QCD 
calculations. 



5.1.3 The integrated dipole functions 



The spin-correlation terms of the splitting functions in Eqs. ( |5.17| ) and ( |5.19| ) have been 

Pj Pij,fi = 0, when = mj,mij = 0. This property en- 



defined such that 



~(m) a 



sures 



11] that the spin correlations vanishes after azimuthal integration. In the evaluation 



of the integrals of the dipole factors over the phase space [dpi{pij,pk)], we can thus replace 
the spin matrices Vj^ fc with their spin averages (Vij^k)- We define 



[(^PiiPij,Pk)] 



{Pi+PjY - ml 



ftg 1 /Anij/ 
2^r(l-e) \W 



(5.22) 



where Uj^ki.^) depends also on the parton masses and Pij ■ pk- In the massless case, the 



integrals Iij,k{^) do not depend on pij ■ pk and become the functions Vjj(e) of Ref. 1TT|. For 
non-vanishing values of the parton masses, these integrals cannot be exactly performed in 
d dimensions in terms of simple functions. We evaluate them by neglecting corrections of 
0(e). Note, however, that we want to compute them for arbitrary (finite or vanishing) 
masses, so we do not want to spoil the commutativity of the massless limit with the limit 
e — s> 0. We thus perform the e expansion uniformly in the parton masses, i.e. the coefficients 
of the 0(e) corrections that we neglect are not singular when one or more parton masses 
vanish. For instance, we do not expand in e terms of the form /ij^*^, whose limits 
and e — > do not commute. 



We decompose the integrals Iij,k{^) into an eikonal part, P^^{e), which contains the soft 
integrals, and a remaining collinear part, J?°)}(e): 



IgQAf^Q^ /^fc; e) = Cf \ 2r^{fiQ, Ilk] e) + I'gQ^ki.l^Q, l^k] e) 



coll 



IQQ,k{^^Q^ f^k;e) — IqQ^ki^^Q, Aifc; e) , 

eik / 



/,,,fc(/ifc; e) = 2Ca 2r^{0, fXk; e) + Igglif^k, e 



(5.23) 
(5.24) 
(5.25) 



The eikonal integral is defined by 
as 1 [A-K^^ 



27rr(i-e)V g2 



{fij,^ik;e) = [dpi{pij,pk)] 



2piPj 1 



5,(1 



yij,k) 



(5.26) 



where mi = and rriij = rrij in the kinematics defined in Sect. p.l.2| . This integral can be 
obtained as a complicated expression of dilogarithmic functions (see App. |^). However, for 
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the calculation of physical cross sections (see Sect. we explicitly need only the symmetric 
part of and that is simpler. We define 



If (/i,-, /ifc; e) = ^ r'\fXj, /ifc; e) ± /"^(/Xfc, /i,-; e) 



reik 



cik / 



SO that 



The symmetric part of the eikonal integral is given by 



J^"^(/ij,/ifc;e) 



where 



-2e 1 



1 - 1,-- 



+ - h - f^k) + 2 Li2(-p) - 2 Li2(l - p) 
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-^Li2(l-p?)-^Li, (l 



Pk 



0(e) 



1 - Vjj^k + 2/iy(l - fij - pI) 
\ 1 + Vij^k + 2/i2/(l - p] - pI) 



{n = j,k), 



'1 - w 



ij,k 



1 + f 



(5.27) 
(5.28) 



(5.29) 



(5.30) 



For the special = 0, which is needed for the g ^ gg splitting, the result for 

I'^^iO, Hk] e) can be easily read off. Using Pj(0, Hk) = 0, Pfc(0, Hk) = Pk and p = 0, we obtain 



Jf (0, /i,; e) = (1 - piV'^, ( 1 - ^-p,"' ) - ^ - ;^Pfc - ^ Li2(l - pD + 0(e) .(5.31) 



2e2 



TT TT 



6 24' 



If both masses vanish, then /i (0, 0) = 0, and we recover the massless result |TT|: 



TT 



/-•^(O, 0; e) = /r (0, 0; e) = ^ - - + 0(e) . 



(5.32) 



For the case of non-vanishing masses, we have pjPk = p, therefore, we can use the relation 

(5.33) 

to obtain the following expansion in e: 



1 - Ipf^ - \pf^ = elnp - e2(lnV, + In^fc) + 0(e=^) 



/+''(Pi,/ifc;e) 



1 



1 



— Inp - Inp In (^1 - {pj + pk) j - 2 Pj ~ 2 ^'^ 

+ ^ + 2 Li2(-p) - 2 Li2(l - P) - ^ Li2(l - p3) - ^ Li2(l - pI) 
+0(e) . 



(5.34) 
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The integrals of the colhnear parts are simpler. The explicit results are 



2e 2e 



2;UQ'^-21n {l-^ik) -^^Q +ln(l-/ifc 



2/ifc(l-2/ifc) 



+ 0(e), 



1 - /ifc I- 11% - III 

-21n(l-/i,)-21n^ ^ 



e ^ ^ V 2 / 3 



2\ 2/ii 



37 l-p^ 



p^ln(^^]-ln^'-^^ 



+ 0(el 



11 50 11 
6e ^ Y ~ y 
+ 0(e: 



1 + /ifc 



.P2 +Pl 



+ ln(l -pfc) 



+ 



1 + Pi ; 1 - p^ 



(2 - 3K)p2 ^ 2pfe 



In 



3(1 -p2) ^1 + ;,, 



with 



Pi 



P2 



r i-pi 



(5.35) 
1 + Pfc 



(5.36) 



(5.37) 



(5.38) 



Finally, we inspect the special case of parametrically-small masses which is relevant 
if one considers massless fermions and collinear singularities regularized by small fermion 
masses, as is often done in electroweak physics. The corresponding asymptotic forms of 
and in this limit read 



rcik 



1 , \ mi 
— , — In — ^ H In — - 

2 



2 ^1 



-In^ 



''ij,k 



ml 



''ij.k 



y ' 



/cou _ i + 1 In ^ 



rcoU 



rcoU 



2 , "^6 
— In— ^ 

3 ^ij,k 

11 50 

67^y ' 



+3 



16 



(5.39) 



where Sij^k = '^PijPk- Note that in these expressions the fermion mass is only regarded as 
collinear regulator. Therefore, on the right-hand side we have neglected not only contribu- 
tions of 0(e), but also all the contributions of 0{rn? / Sij^kji even in coefficients of e poles. 
The collinear logarithmic contributions in the Q gQ splitting are related to those in the 
fermion-photon splitting considered in Ref. |]25|, where soft divergences were regularized 
by an infinitesimal photon mass. 
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5.2 Final-state emitter and initial-state spectator 



The dipole contribution V^ - (see Fig. |1|) to the factorization formula ( |5.1| ) is defined in close 
analogy to that in Eq. (|5.2j ): 



ij,- ■ ■;a,- ■ -l "2 • • • ,U, • • • ; a, • • ■)m,a ■ (5.40) 



The m-parton matrix elements on the right-hand side are obtained from the original (m+1)- 
parton matrix element by replacing the final-state partons i and j with the emitter parton 
ij and the initial-state parton a with the spectator parton a. 

The final-state parton momenta pi and pj have arbitrary masses, the initial-state mo- 
mentum Pa is massless and their total transferred momentum is denoted by Q (not to be 
confused with Q in Eq. (|5.3|) ), 



Pt =T^i, Pj =^j^ Pa = 0, Q=Pi+ Pj - Pa- (5.41) 

The on-shell mass of the emitter is rriij, with rriij < rrii + rrij. 



As discussed in Sect. |2.2| , we do not consider massive partons in the initial state in this 
paper. A treatment of massive initial-state particles is described in Ref. for photon 
emission off fermions. 



5.2.1 Kinematics and phase-space factorization 



We introduce the auxiliary variables 



PaPi + PaPj - PiPj + 



2 



PaPi + PaPj 



PaPi 



^0 



PaPj 



PaPi + PaPj PaPi + PaPj 

and we define the momenta of the emitter and spectator as 

Pa = a;y,aP^, P% =Pi + Pj - (1 - Xij,a)p'^. 

These momenta obey their mass-shell relations and momentum conservation, 

pI = Pa = 0, Q= Pij 

We also define the rescaled parton masses /i„ as 



Pa- 



rrir. 



2pijPa 



(5.42) 



(5.43) 



(5.44) 



(5.45) 



The momentum definition in Eq. ( ^.43|) coincides with that of Ref. for the specific case 
(with rriij = f^j and rrii = 0) considered there and with Ref. |^ for rriij = nT'i and rrij = 0, 
but differs for arbitrary combinations of masses. 
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Since the introduction of the rescaled momentum pa = Xij^aPa defines a new (boosted) 
frame, the factorization of the two-particle phase space d(j){pij]Q + pa) from the three- 
particle phase space d(f){pi,pj] Q + Pa) involves a convolution over a boost parameter x, 

d(p{pi,pj; Q + Pa) = [ dx d(f)(pij; Q + xpa) [dpi{pij;pa, x)] 9(x+ - x) , (5.46) 



with the upper limit 



(5.47) 



The single-parton phase space [dpi{pij;pa,x)] is 

1 

[dpi{pij;pa,x)] = -(27r)"^+^'(2pijPa)^~'y^ dxij^aS{x - Xij^a) (1 -x + i^^y 



z+ (x) 



X / d'^-'n I dz^ [z+ix) - ~ZiV [zi - 



with the limits 



1 - a; + /i-j + /if - /i| ± ^(1 - X + - fif - - 4/if/i2 

2{l-x + fi%) 



(5.48) 



(5.49) 



The d'^ integration is over the solid angle Q perpendicular to pij and pa- 



5.2.2 The dipole splitting functions 



We give the functions V"^- in Eq. ( p. 40 ) for the two different QCD splitting processes 



Q giPi) + QiPj)- mi = Q and rrij = rriij = mg, 
g Qipi) + Q{Pj)- rrii = ruj = niQ and niij = 0. 



The case Q —>■ gQ is formally identical to Q ^ gQ- The splitting g gg is already given 
in Eq. (5.40) of Ref. 0], since no massive parton is involved in this case. Specifically, we 
define the dipole functions 



z ^ ~ Zj 

^ Xij a Zj 



m 



Q 



Pipj 



e(l - Zj) \ 6ss 



(/^|V^q|z/) =8V^«sTr^-(7 



ZiPi - Zjpfj ZiPi - ZjPj 



(5.50) 
(5.51) 

(5.52) 



(Pi+Pj)^ 

(^qq) = 8vr/i2^asTR|l - y^(^+ - z.^ - zJ) ^ 
with z± given in Eq. (|5.49|) . 

In App. |C|, we give the dipole splitting functions that are relevant for SUSY QCD 
calculations. 
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5.2.3 The integrated dipole functions 



The integral of the spin- averaged dipole functions (V"^) over the singular phase space 
[dpi{pij;pa,x)] is defined by 



J [dp i{Pij;Pa,x)] 



(Pi + Pj] 



in) 



27rr 1 



'^PijPa 



(5.53) 



where /f • (x; e) depends also on pij ■ Pa and the parton masses. In the massless case, /f • does 



not depend on pij ■ pa and becomes the function Vjj(x; e) of Ref. ITT 



For e = the function Igqix; e) is well-defined everywhere, except at the endpoint x = 1 

where the integral is singular as ^^^jz^- To avoid this singularity, the limit e — has to 
be performed uniformly in x, which is achieved by decomposing Igqix; e) into a continuum 
and an endpoint part using the usual '+ '-distribution, which is defined by 



dx (f{x)) g{x) = f dx f{x) [g{x) - g{l)] . 



, (5.54) 

In this manner, Igqix; e) defines an x-distribution whose coefficients contain poles in e. 

In the case of the splitting process g QQ, the endpoint is located at x+ = 1 — 4/1^, 
where the function lQQ[x]e) is not singular. However, it becomes singular as y~: for 
vanishing quark masses. To ensure a smooth massless behaviour, the limit e has to be 
performed uniformly both in x and in the quark mass. To achieve this, we have to use a 
modified version of the '-|- '-distribution, which we call '-distribution, 



dx (/(a;))_ 



dx /(x)0(x+ - x) [g{x) - g{x+)] . 



(5.55) 



The 'x+ '-distribution reduces to the usual '-|- '-distribution by setting x+ = 1 and, in par- 
ticular, it tends to the '-|- '-distribution smoothly when /ig — >• in Iqq{x; e). 

After decomposing I°j into a continuum and an endpoint part, we further decompose the 

endpoint part into a singular ( Jj"'^) and a finite (Jij^^) piece, such that the latter vanishes 
for vanishing masses: 

PgQix; e) = Cf {[J^q{x, /xq)]+ + 5(1 - x) [j;|(/iQ; e) + J^'^^'ifiQ)] } + 0(e) , (5.56) 
PQQix; e) = Tr {[J^q(x, /xq)],.^ + 5(x+ - x) [j5|(/iQ; e) + r^fif^Q)] } + 0(e) ,(5.57) 

where x+ = 1 — 4/iQ and Jij^^ifJ^q = 0) = . 



The three contributions to Eq. ( |5.56| ) are 

/ 1 — X 



1 — X 



2(1- x + Illy 

+ (^] ln(2 + /i^-x 



l + ln(l-x + /iy 



.1 — X 



TT 

"6 



TT 



-ln(l + /i^) + - + ---, 



(5.58) 
(5.59) 
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(5.60) 



Note that the notation in Eq. (|5.58|) is symbohc: [JgQ^x, fiQ)]+ contains all the plus dis- 
tribution contributions, but is not itself a plus distribution. The same is true of [Jgg{x)]+ 
defined in Eq. ( |5.67| ). 

Setting the mass of the quark to zero, we immediately obtain the corresponding expres- 
sion in the massless theory (see Eq. (5.57) in Ref. [jlll). If the quark mass is non- vanishing, 
we can expand the singular part of the endpoint contribution as 



1,92 ^ "^1 ^ ^2 

m Un H \ — In u7) H tt 



(5.61) 



The three contributions to Eq. (|5.57| ) are 



x+ 



2 / 1 - X + 2^^ I 4/i| 



(1 -x)2 
2 



3e 



1 - ^^Q' 



1 — X 

10 



10 



1 + ^/l - 4^ 



(5.62) 

(5.63) 
(5.64) 



The massless result of Eq. (5.58) in Ref. [|TT| is again recovered by simply setting mq = 0. 
If mq 7^ 0, the coUinear pole 1/e is replaced by the logarithm of the scaled quark mass. 



> 0) = -gln/^Q - Y • (5.65) 
As mentioned above, the case of the splitting process g ^ gg is already discussed in 

(5.66) 



Ref. [|TT]. We recall the result by using the present notation: 



2Ca {[J;,(x)]+ + 5(1 - x) J,f (e)} + 0(e) 



where 



2 In. ' 



11 1 



1 — X 1 — X 61 — X 

1 11 67 TT^ 

6^^ 18 ~y ■ 



1 — X 



ln(2 - x) 



(5.67) 
(5.68) 



Finally, we consider the asymptotics in the case of a parametrically-small quark mass 



niq = nq J2pijPa, uscd as coUinear regulator: 



ln(l — x) 



2 VI - X 



+ 



+S{l-x)(- (in/ij + l) - i InVg + ^ Infil - ^tt^ + 



1 — X 
1 
2 



ln(2 - x) 



(5.69) 
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T, 



R 



X 



5(1 



2,2 10 



(5.70) 



Analogously to Eq. ( ^.391) , on the right-hand side we have neglected all the contributions 
of 0(/iQ), independently of their e dependence. Comparing these results with Eqs. (5.57) 
and (5.58) in Ref. |]ll|, we notice that the continuum parts reduce to those of the exactly 
massless expected since all the singularities are related to the final-state emitter. 

The 1/e coUinear poles of the massless theory in the endpoint parts are replaced by the 
mass logarithms, which in the Q gQ splitting are related to those of the fermion-photon 
splitting p5[1 . 



We conclude this section with a comment on a subtle point related to the definition of 
the x-distributions or, more precisely, to the definition of the space of test functions onto 
which the x-distributions act. Indeed, the integrals Ifj depend on two variables, x and pij, 
which are both integration variables in the calculation of the cross section. Our definition 



in Eqs. ( ^.56|) and ( ^.57|) , in terms of continuum and endpoint parts, considers /f- as an 



x-distribution that acts on functions of the variables x and Pij, and Pij is regarded as an 
additional and independent integration variable that has to be kept fixed in the subtraction 
prescriptions of Eqs. ( |5.54| ) and ( |5.55| ). However, the x-distributions can also be defined in 
a different way, for instance, by regarding Q = pij — xpa (instead of pij) as the independent 
integration variable. This alternative definition is explicitly worked out in App. where 
we show that the replacement pij —^Q= pij — xpa, being x dependent, has to be handled 
with care in the definition of the x-distributions. 



5.3 Initial-state emitter and final-state spectator 

We define the dipole contribution (see Fig. |^) to the factorization formula ( [5.1| ) in close 
analogy to that in Eq. (|5.2| ): 



V]\pi,...,Pm+i;Pa,---) = 

^ ^ .(..., j,. .. ..|^47^^Vr|... J,. ..;^,...L-.. (5.71) 



9n n- 'T . • rn,ai^' • • ' J ' • • • ' ) • • • I rp2 j 

The m-parton matrix elements on the right-hand side are obtained from the original {m+1)- 
parton matrix elements by replacing the initial-state parton a with the emitter parton ai, 
the final-state parton j with j and dropping the final-state parton i. 

The final-state parton momentum pi and the initial-state momentum pa are massless 
and the mass of spectator momentum pj is arbitrary. The total transferred momentum is 
denoted by Q (not to be confused with Q in Eq. (|5.ij|)), 

p1 = 0, p] = m], Pa = 0> Q = Pi + Pj - Pa- (5.72) 

The on-shell mass of the emitter is m„i = 0. 
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5.3.1 Kinematics and phase-space factorization 



We can completely take over the kinematics of Sect. |5.2.1| after replacing pa Pai, Pij Pj 



rrij — >■ and rriij — >■ rrij there. Thus, the new momenta read 

Pai ~ ^ij,aPa^ Pj=Pi~'t~Pj^i^~ ^ij,a)Pai (5.73) 

where 

^ ^ PaPi + PaPj - PiPj ^ _ PaPi ^ _ PaPj 

^^'^ PaPi + PaPj ' * PaPi + PaPj ' ^ PaPi + PaPj ' 

They obey the mass-shell relations 

pI=pI= Pli = 0, p^j = p] = m] (5.75) 
and momentum conservation 

Q = Pi + Pj - Pa = Pj - Pai- (5.76) 

The rescaled mass fij of the parton j is defined by 

/z, = . (5.77) 

'2pjPa 

For the sake of completeness we also write down the phase-space factorization, 

J d(f){pi,pj]Q + Pa) = dx J d(f)[pj;Q + xpa) J[dpi{pj]Pa,x)] (5.78) 

with 

[dpi{pj;pa,x)] = -{2ny^^^'{2pjPay~^ / dxij^aS{x - Xij^a) (1 - a; + /i^)~' 
4 JO 

X J d'^-^n ^'^^''^ dzi [z+{x) - Zi]" Z-' , (5.79) 
where the upper limit of the 2;-integration is 

1 — X 
1 — X + /ij 



z^{x) = (5.80) 



5.3.2 The dipole functions 

We consider the splittings q —>■ qg ot gq, g —>■ qq or qq, g gg of massless partons in the 
presence of the massive spectator parton j. The corresponding dipole functions read as 
follows: 

{s\yf\s') = STTfi^'a^Cp I — ^ _ , - 1 - Xij^a - e(l - Xij^a) \ Sss' 

= iyf)5ss' , (5.81) 
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{s\\f\s') = 

(/.|V|V) = 
(Vf) = 



(Vf) 



(Vf , 



1 - Xij^a 'iZi^j fp'^ p^j \ fp'^ p^j 



Xij^a PiPj ^ Zi Zj J \ Zi Zj 
2(1 — Xij^a) 



1 



2 Xij^a Zj 



1 ~l~ 2Jjj^a(l Xij^a] 







(pl 








\Zi 





167r//^^Q;sCA 



2 



1 ~l~ •^ij,a(l •^ij,a) ~l~ 



(5.82) 
(5.83) 

(5.84) 



(5.85) 



2 ~ 



1 Xij,a f^i Z'l 



■^ij,a Zj 



(5.86) 



5.3.3 The integrated dipole functions 

We define the integral of (Vf) over the singular phase space [dpi{pj;pa,x)] by 

mpr,Pa,x)] ^ ^ (Vf ) ^ ^^,,^(^Ylf{x,^i,-e), (5.87) 



2paPi Ws(a) 



27rr(l-e) V2p,p„ 



where the factor ns(ai)/ns{a) turns the spin-average over ai into the one over a. Splitting 
off the massless part Ij^{x, 0; e), the results for /"^(x, fij; e) can be written in the compact 
form 



If{x, ,ij- e) = If{x, 0; 6) + P,tg(a;) In - 



1 -X 



X + fjL^j 



.2 2/^1 



+ d,bT'i ^ In ( ^-^ ^ 1 + 25'^^ ( ^— 



In 



2-x 



+ 5"^5(l-a;) 
+ 0(e), 



- ln(l + /i •) + - ln2(l + /x^) + 2 Li2 



l + li 



(5.88) 



where the auxiliary functions P^egl^) given by 



PZi^) ^ -C^{1 + x) , 



l + (l-x)^ 



P,%{x)^T^[x' + {l-xf 
l-x 



X 



Pilix) = 2Ca 



X 



— 1 + x{l — x) 



(5.89) 



Since the kinematics of the spectator parton j is not connected with any singularity, the 
limit rrij — > is smooth in all cases. The massless part can be written in the following 
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compact form, 



/f(x,0;e) 



P''\x) + P^fjx) ln(l ~x) + 2 6'"' Tl 



ab /T-i2 



'ln(l — a;) 
\ — X 



ln(2 - x) 

\ — X 



(5.90) 



where 7a are the following flavour constants: 



7. = 



(5.91) 



and P'"^(x) is the e-dependent part of the d-dimensional Altarelli-Parisi splitting functions 
in Eqs. (TO - CT) : 



{Pab{x; e)) = {Pab{x; e = 0)) - e P"'\x) + 0{^) . (5.92) 

Explicitly, 

P"^''{x) = P'^^{l-x) =Cf{1-x) , P'^%x) =2T^x{l-x) , P'55(x) = 0. (5.93) 



The functions P°'^{x) in Eq. (|5.9CI|) are the usual regularized Altarelli-Parisi functions that 
can be written as [ITTl 



P'^'ix) = P,%{x) + 5' 



ab 



(5.94) 



5.4 Initial-state emitter and initial-state spectator 



Since we allow for massive partons only in the final state in this paper, the case of emitter 
and spectator both from the initial state is already completely covered in Ref. [^, where 
all cases for massless QCD partons are presented. Results for the Q Qg splitting with 
a massive quark in the initial state can be deduced from the treatment of photon emission 



off fermions described in Ref. [25| 



6 QCD cross sections at NLO 



In this section we describe our method for evaluating QCD cross sections of massive and/or 
massless partons. We start with the precise definitions of the three terms on the right-hand 
side of Eq. ( |2.13| ), which can be considered the final expressions, ready for implementing 
in a partonic Monte Carlo program. Then we derive the explicit forms of the insertion 
operators I, P and K relevant to the three different cases of initial state: elementary 
particle reactions involving zero, one or two initial-state hadrons. More details on our 



notation can be found in Sects. 6-11 of Ref. 11 . 
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6.1 Final formulae 



In order to present the final formulae of our algorithm, we start with the most general 
case of two incident hadrons and then point out the simplifications if one or both of the 
incoming particles are leptons. 

In the case of processes with two initial-state hadrons carrying momenta and p^, 
the calculation of the QCD cross section requires the introduction of parton distributions. 
If we denote by fa/AijIi f^p) the density of partons of type a in the incoming hadron A, the 
hadronic cross section is given by 

(^{PA,Pb) =Yj f ^Va fa/A{Va, fJ^p) f dr]b fb/B{Vb, I^'f) CTabiVaPA, VbPB) , (6.1) 
a,b ^0 ^0 

where aab is the partonic cross section, which in the first two orders of the perturbative 
expansion has the schematic form given in Eq. ( |2.1|) , 



(TabiPa^Pb) = (Tab iPa,Pb) + (Tab i.Pa, Pb] f) ■ (6-2) 

If one or two of the colliding particles is pointlike, i.e. a lepton in the Standard Model, 
then the integration becomes trivial by replacing the corresponding parton density by 
5(1 — i])6aA- The LO parton-level cross section is given by 



(Tal^{Pa,Pb) = / da^^{pa,Pb) 
J m 



d^^^\pa,Pb) . , \Mm,ab{{P^,m}■,Pa,Pb)\'F!^\{p,}■,Pa,Pb) ,(6.3) 

where a and b denote the flavours of the incoming partons, Pa and pb are their momenta 
and nc{a),nc{b) are their number of colours (for an incoming lepton Hc is set to 1). The 
matrix element |A1m,afeP is the square of the tree-level amplitude to produce m final-state 
partons. The factor 

1 1 

d$('")(p„Pfe) =Afin -J2d(Pm{pi,...,Pm;Pa+Pb)^ (6.4) 

contains (i) the averaging over the number of initial-state polarizations {ns{i) number of 
states for parton i, i = a or b), (ii) the flux factor J-'{pa ■ Pb), (hi) a summation over all 
configurations with m partons in the final state, (iv) the Lorentz invariant phase space of 
m final-state partons d(f)m, (v) a Bose symmetry factor S{m} for identical partons in the 
final state and (vi) A/^n, which includes all factors that are QCD independent. Finally, the 
function Fj^'' in Eq. ( |6.3| ) defines the jet observable we want to compute. It has to fulfil 
the formal requirements of infrared and collinear safety. Furthermore, in order for our cross 
section to have a smooth small-mass limit, Fj should also be quasi- collinear safe. 

The general definition of quasi-coUinear safety is analogous to that of collinear safety: 
if a pair of partons i and j whose masses and relative transverse momenta are all 0{XQ), 
where Q is the hard scale, then as A — 0, Fj^'^^\. . . ,Pi,Pj, . . .) must tend to the value of 
Fj^\. . . ,pi + pj, . . .). Any collinear-safe observable for massless partons can be extended 
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to be quasi-collinear-safe for massive partons. However, there is a subtlety related to the 
flavour-dependence of Fj. To respect infrared and coUinear safety, an observable must not 
depend on the identities of massless partons in the final state. It can however depend on 
the identities of massive partons, for example the thrust distribution in e^e~ annihilation 
events that contain at least one bottom quark is an infrared- and collinear-safe quantity. 
However it is not quasi-collinear safe, because events in which the hh pair is quasi-collinear 
contribute but identical events with the hh pair replaced by a gluon do not. Careful defini- 
tion of the observable can render it quasi-collinear safe and hence finite in the high-energy 
limit. One solution in this example would be to only include events in which the h and 
h lie in opposite thrust hemispheres. Our method is perfectly capable of calculating cross 
sections for observables that are not quasi-collinear safe, but our careful treatment of the 
quasi-collinear limit only guarantees improved numerical convergence for quasi-collinear 
safe observables. 



According to the general notation in Eq. (|2.13|) , the full NLO contribution is a sum of 
three terms, 

(^a6^°(Pa,P6;/i|) = (^Ib^^'^^^^iPa^Pb) + (^T^ ^"^^ {Pa, Pb) 

, f^, r^NLO{m}/ 2\ , -NLO{m}/ 2\] (ar:\ 

+ dx (T^fc {x;XPa,Pb,^iF) + (^ab {x] Pa, XPb, ^p) ■ (6-5) 

Jo L J 



The first contribution has (m + l)-parton kinematics and is given by the following 
expression 



NLO \m+l} 
(^ab 



{Pa,Pb) 



m+1 



(d0-S,(Pa,Pb) 



e=0 



i dipoles 



€=0 



Jd^^"'+'\pa,Pb) 



nc{a)nc{h) 



\Mm+l,ab{{Pi,^i}',PaiPb)\'^ Flf'~^^\{Pi};Pa,Pb) 



Y (l^-F(-))({p„m,};p„p,) 

dipoles 



(6.6) 



where Airn+i,ab is the tree-level matrix element with m + 1 partons in the final state and 
I^dipoies ■ -F'-™'') {{pi, mi}]pa,Pb) is the following sum of the dipole functions: 



J2 J2 Vij,k{pi, . . . ,Pm+l]Pa,Pb) Ff'\pi, . . .pij,Pk, . . . ,Pm+l]Pa,Pb) 

pairs k^ij 

+ Y [^ij (Pl' • • • ' Pm+i;Pa, Pb) F^r^ {pi,...Pij,..., Pm+l]Pa, Pb) + {a ^ h) 

pairs 

+ [^kiPU - ■ ■ ,Pm+i;Pa,Pb) F^'^\pi,. . .pk,. . . ,Pm+i;Pai,Pb) + (fl ^ 6) 

i ky^i 

+ J2 [V''''\pi,...,Pm+i;Pa,Pb) F^"'\pi,. . .,p^+i;Pai,Pb) + (o ^ 6) 



(6.7) 



where (a ^ h) indicates that the roles of a and h are reversed. For example, in the 
second hne, Vf^ becomes V^j and Pa,Pb becomes Pa,Pb- In the last line of Eq. ( |6.7| ), the 
function Fj^^ depends on the final-state momenta pi, . . . , Pm+i, which are obtained from the 
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original momenta Pi, ■ ■ ■ ,Pm+i by a Lorentz tranformation (see Sect. 5.5 and the Erratum 
in Ref. ITl). 



If an incoming particle is a lepton, then the terms with I) carrying the corresponding 
parton-flavour index are absent from the sum, i.e. for one initial-state hadron all terms 
corresponding to the (a <-h> b) interchange and the full last line has to be dropped and for 
zero initial-state hadrons, only the first line remains. 

The NLO contribution with m-parton kinematics is obtained by adding the virtual cross 
section and the singular part (with some accompanying finite terms) of the integral of the 
subtraction term over the factorized one-particle phase space. The former is given by the 
interference term between the tree and one- loop amphtudes | A^m,at)|(i_ioop) latter 
in terms of an insertion operator /(e): 



^NLO {m} ^^^^ ^ / r^l^V (p^^ + da^hipa, Pb) ® I 
J m 



e=0 



+m,a6(l, . . . , m; a, 6 I J(e) |1, . . . , m; a, h)m,ab 



F!r\{p.};Pa,Pb),m 



t=0 



where the operator /(e) is spelled out explicitly for the various cases of interest: in 
Eq. ( |6.16| ) to be inserted into ^(1, • • • , "^1 11, • • • , '^)m (case of no initial-state hadrons), 



( |6.51| ) into m,a(l, . . . , m; a| |1, . . . , m; a)m,a (case of one initial-state hadron) and ( |6.66|) into 



m,ab{^, . . . , m; a, 6| |1, . . . , m; a, b)m,ab (case of two initial-state hadrons). 

The third term on the right-hand side of Eq. ( |6.5D is a finite remainder left after can- 
cellation of the e-poles of the collinear counterterm. It is a sum of as many contributions 
as the number of initial-state hadrons, each obtained by integrating a cross section with 
m-parton kinematics with respect to the fraction x of the longitudinal momentum carried 
by one of the incoming partons. When this parton is the parton a, we explicitly have: 

/'da;a,T^'"^(x;a;p„p,,/i|) = 5:/'da; / [da,\(xp„p,) ® (/C + P)"'»'(x)] 

= ^ f'dx /d$('")(xp„p,)Ff)(te};xp„p,) 

x^yb(l, . . . , m; xpa,Pb\ /^"'"'(x) + P"''''(a;,/i^) |1, . . . ,m; xpa,Pb)m,a'b , (6.9) 
where the colour-charge operators P and K are respectively defined in Eqs. ( |6.53| ) and 



( |6.55| ) to be inserted into m,a' ■ ■ ■ a^Pal |1, ■ ■ ■ ,m; xpa)m,a' and Eqs. (|6.67|) and (|6.6S 



into m,a'b(l, • • • ■,TTT']xpa-,Pb\\^i ■ ■ ■ i ^] ^Pai Pb) m,a'b- If there are two initial-state hadrons in 
the collision, then we also have to add the integral of a^l''^^^^ {x-^pa-,xpbi ^ip)-, which has a 
completely analogous expression to Eq. (|6.9[), apart from the replacements xpa — > PaiPb 
xpb and J2a' J2b'- 



Equations (|6.3|) , ( |6.6|) and ( |6.9| ) are directly evaluated in four space-time dimensions. 
As for Eq. ( |6.8|) , one should first cancel analytically the e poles of the one-loop matrix 
element with those of the insertion operator /, perform the limit e — and then carry out 
the phase-space integration in four space-time dimensions. 



35 



In the following subsections we derive the explicit form of the /, K and P inser- 
tion operators, for the various specific cases, relevant to lepton-lepton, lepton-hadron and 
hadron-hadron scattering. 



6.2 Jet cross sections with no initial-state hadrons 



In processes with no initial-state hadrons, the QCD cross section for jet observables at 
NLO accuracy is given by Eqs. and (|2.2| ), and the NLO contribution is rewritten 

into the form of Eq. ( p.3|) . In order to calculate the integral of the subtraction term over 
the factorized one-parton phase space, we start by writing J^^i da^ using the notation 
introduced in Eqs. ( [OD , (|6.4| ) and the phase-space factorization of Eq. (|5.10|) . We find 
that the integral of the subtraction term can be written in the following form: 

/ da^=Xin J2 J2 J2 d(l)m{Pl,---,Pij,---,Pk,---,Pm+i;Q) 

J m+1 r ,11 i.-l; ; Jnri 



{m+1} pairs kj^i,j 

i,3 



X-^r^ F'f\pi, ■ ■ ■,Pij, ■ ■ ■,Pk, ■ ■ • ,Pm+l) 

•^{m+l} 

X [dpi{pij,pk)]Vij^k{pi, ■ ■ ■ ,Pm+i) ■ (6.10) 



For the integration over dpi we write the dipole factor as in Eq. ( p.2| 
T^ijAPi^ ■ ■ -^Pm+i) = 



Vij^k 1 



— TPT l-^m^(Pl> • • • ^Pij^ . . . . . . ,Pm+l] 

■(6.11) 

where Ai^J^'^ is a colour-correlated m-parton amplitude (see Eq. ( |3.4| )) depending only on 
pi, . . .pij,Pk, ■ ■ • ,Pm+i while Vij^k/iiPi + Pjf - mfj] depends only on Pi,Pj,Pk or, equiva- 
lently, Pi,Pij,Pk- The subscript h on the square bracket in Eq. ( |6.11| ) means that Vij^k and 
are still coupled in helicity space. According to the discussion in Sect. |5.1.3|, this 



helicity-space coupling vanishes after integration over [dpi{pij,pk)], and the integral of the 
subtraction term remains the same if we exchange the splitting matrices for the azimuthally 
averaged splitting matrices. Then we can perform the counting of the symmetry factors in 
order to change S^m+i} to S^m}, which follows exactly the procedure of the massless case[|T^] 
leading to Eq. ( |6.10|) rewritten in a completely factorized form of an m-parton contribution 



multiplied by a singular integral: 

/ do-^ = -Afin I #m(Pl, ...,Pm]Q) ^j"'Hpl, ■ ■ ■ , Pm) 

Jm+1 Jm O^m} 

X E E \ML\P^, ■ ■ ■ ,Pm)P (/ [dpM,,Pk)] jjy^ , ) . . (6.12) 

k^k 

The substitutions in the last factor are to be performed after the integral is calculated 
using the results in Sect. p.l.3| . In this context we observe that the colour-connected m- 



parton amplitude, A4{n symmetric in its upper indices, and we sum over all pairs of j,k 
in Eq. ( |6.12| ). Consequently, the antisymmetric part of the eikonal integral, I'^^ does not 
contribute to the integral of do"^. 
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After these manipulations, the final result for /^^^dcr^ in Eq. ( 6.12|) can be written as 
follows, 



m+l 



dcr^ (g) Im{e,^^; {pi,mi}) 



(6.13) 



where the notation 



on the right-hand side means that one has to write down 
the expression for dcr^ and then replace the corresponding matrix element squared at the 
Born level 



\M 



ml — m) 



11, . . . ,m||l, . . . ,m), 



by 



^(1, . . . ,m| /„^(e,/i^; {pi,mi}) |1, . . . ,m)^ , 



(6.14) 



(6.15) 



where the insertion operator 1^ depends on the colour charges, momenta and masses of 
the m final-state partons in da^: 



.2 / 



TC 



Vj{sjk, rrij, TUk, {mp}; e, k) - — + Tj{fi, ruj, {mp}; e) 



u 

+7j In + + Kj + 0(e) 



(6.16) 



where Sjk = 2pjPk and {nip} denotes the set of masses of those massive quarks that may 
appear in the g ^ qq splitting. The constants 7(j, were defined in Eq. ( p.91|) , while Ka 



denote the constants introduced in Ref. 

/7 7^2^ 



6 



18 6 



9 



(6.17) 



In Eq. ( |6.16| ) the Vj kernels depend on the fiavour of parton j and on the momenta and 
masses of both partons j and k. Thus colour correlations are accompanied by complicated 
mass correlations. We decompose these functions into a sum of two contributions: the 
first is symmetric with respect to interchange of the indices j and k and singular in four 
dimensions or for vanishing masses and the second is neither symmetric nor singular, 

Vj{sjk,mj,mk, {mp}; e, k) = V^^\sjk,mj,mk; e) + Vf^\sjk,mj,mk, {mp}; k) . (6.18) 

The singular terms have the following form: 



1 [Qjk 



V^^\sjk,mj,mk; e) 



(6.19) 

where we introduced the abbreviation = Sjk + + m\ and Q{x) is the step function 
with 6(0) = 0. Here and in the following the quantities p, pj and pk are obtained from 
Eq. (|5.30|) with the substitutions p\ m,^/Q|^ (n = j, k) and Vij^k Vjk- In practice the 
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expanded forms of V*^^^ are needed, 

V^^\sjk,mj > 0,mfc > 0;e) 



1 

Vjk 



1 1 2 2 "^122 

Inp - - In pj - - In 



4 



77" 



H In p In 



1 



1 , 1 , , vr^ 
— H In ^ In^ ^ 

2e2 2e Sjk 4 12 



-- In — - In ^'^^ 



Sjk 

llniln ^^-^ 



2 Sjk 2 Qjjt Qj-fc 



2 ' 



.2 ■ 



(6.20) 



Of course, the decomposition defined in Eq. ( |6.18D is not unique. One can still shift finite 
terms between the non-singular and singular terms. In particular, in Ref. |^ we have not 
included those terms in the singular parts that vanish for vanishing masses (terms in the 
second lines of the equations in ( |6.2(]| )). 

The non-singular terms depend on the flavours and masses of partons j and k. If j is 
a quark, then V^^^^ does not depend on {mp} and so we suppress these arguments. If 
both j and k are massive quarks, then 



Sjk,mj,mk) 



Iq Sjk 



In 



+ — 

Vjk 

+ In 



lnp2 ln(l + p') + 2 Li2(p') - Li2(l - p]) - U^il - pi 



~6 



Qjk - ruk 

Qjk 

rrik 



2 In 



{Qjk - rukY 



1m] , 
m 



jk 



Sjk Qjk — "^fc 



Qjk — ^k 



^ 2mfc(2mfc - Qjk) ^ vr^ 



Sjk 



(6.21) 



where we used the abbreviation Qjk 
parton, this reduces to 



Qj^. If i is a massive quark, but /c is a massless 



Vr\s,k.m,^^) = ^,\r.^ + \-U, 



Q 



Sjk_ 



_2 1n^-:^ln 



q -^jk ^ \^jk/ 

If j is a massless quark and A; is a massive parton, then 

Qjk - rUk 2mk 



Sjk 



Q' 



Vfs)(.,„0,m,) = ^ 



In 



Sjk 



Q' 



21n- 



jk 



Qjk 



Qjk + mk 



jk 



^ " Til 



(6.22) 



(6.23) 



If j is a gluon, then 

V'j^^\sjk,0,mk, {mp}; k) 



73 



In 



Sjk 



Q' 



2 In 



jk 



Qjk — ^k 
Qjk 



Qjk + rUk 



+ 



Li2 



Sjk 



jk, 
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F=l 

+ -— 2^ In 



3Ca 



1 Qjk — ruk , 
In „ h 



Qjk + mk 



+ In 



1 +Pi 



Pi 



(3 + P?) 



1 , m 
- m 
2 



ml 



3Ca 



Q 



2 

aux 



+ K 



2\ ml 



3/ SjJ; 



1 In 



2mi 



F=l 



P2- Pl 
,P2 + Pl 



-In 



i-pA i 

1 + pJ 



(6.24) 



where Np is the number of heavy flavours and Np' is the number of those heavy flavours 
for which Sjk > AraF{raF + mk)- Note that for mi;' — 0, we can set Np^ = Np because the 
region of vanishing sjk gives a vanishing contribution to infrared-safe observables. Thus 



,mk,{mF};K) is finite in the massless limit. The auxiliary mass scale Qaux has 
been introduced only to ensure this finiteness property. It can be chosen arbitrarily|U] and 
will cancel against a similar contribution in Tg{{mF}; e) (cf. Eq. ( |6.27| )). In Eq. ( |6.24] ) pi 
and p2 as defined in Eq. ( ^.381) can be rewritten as 



Pl 



Arrip 



\ {Qjk -rUkY 



P2 



Amp 



\ Q]k - < 



(6.25) 



We see that choosing n = 2/?) simplifies Eq. ( |6.24|) considerably. 

Finally, if both j and k are massless partons, then V^'^^''(sjfc, 0, 0) = and Eq. ( |6.24| ) 
reduces to 



Vfs)(s.fc,0,0,{m^}) 



N- 



■jk 



3C ^ 



A F=l 



In 



1 + Pi 



f(3 + p;)-liu!!^ 



2 
F 

Sjk J 



2Tr 

3Cl 



'^A ^1 Q 



mp 



2 

aux 



(6.26) 



where we dropped n from the arguments because the dependence on n in Eq. (|6.24| ) is 
proportional to the mass of parton k. If the heavy partons are completely absent, then 
V^^^\sjk, 0, 0, {}) = 0, with {} meaning the empty set. 

The functions Tj depend on the flavour of the parton j and on the parton masses. In 
the case of gluons and massless quarks (antiquarks) we have (with dummy arguments 
suppressed) 



1 2„ ^. m^, 



Tg{{mF}; e) = 7 73 - 3 II In ^ 



r.(6) 



2 

aux 



7g 



(6.27) 
(6.28) 



while for massive quarks (antiquarks) we find 

.2 \ 



m 



m 



r,(/i,m,;e) = - - In ^ - 2 + 7, In ^ = 



P 



"in Ref. pi[ Qaux was chosen to be the renorniahzation scale. 



^In^ 

2 /i2 



(6.29) 
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In Eq. ( 6.27] ) the summation runs over the heavy flavours. 



Using formulae (|6.16|) -( |6.29|) , it is easy to verify that the /^(e, /x^; {pi,mi}) operator 
tends to the corresponding operator of the massless theory smoothly, 



lim /^(e,/i^; {pi,mi}) = /(e) + 0(e) 

{mi}-*0 



where /(e) is defined in Eq. (7.26) of Ref. 



(6.30) 



The NLO QCD cross section in Eq. ( p.3|) is finite for infrared-safe jet observables. We 
constructed the subtraction terms such that the first term is not only finite independently of 
the parton masses, but it also tends to the corresponding massless expressions smoothly if 
Fj is quasi-collinear safe (see Eq. (|2.6|) ). Therefore, the second term, spelled out in Eq. (|6.8|) , 
also has to be finite and, using Eq. ( |6.30| ), tends smoothly to the corresponding massless 
expressions in the zero-mass limit. Thus all the e poles and In(mj) terms that are singular 
in the zero-mass limit in do"^ must be cancelled by corresponding terms in do"^ ® Im{^), 
which enables us to write down the singular terms of the one-loop amplitudes, as spelled 



out explicitly in Ref. [31 . 



6.3 Jet cross sections with one initial-state hadron 

As mentioned before, in this paper we consider only massless partons in the initial state, 
therefore, the implementation of the subtraction procedure for processes with hadrons in 
the initial state is the same as described in Ref. |]TT| (see also Sect. (|6.1|)). The presence of 
the parton masses induces non-trivial changes in the evaluation of the sum of the integral 
of the full subtraction term and the collinear counterterm, 

/ dcr^iPa) + I da,C(p,) , (6.31) 

where the index a refers to the flavour of the incoming parton and Pa to its momentum. 
The full subtraction term da^{pa) is a sum of three terms, corresponding to the first three 
terms in Eq. (|5.1|) , 

da^ip^) = da^'ip^) + dafip^) + daf{p^) . (6.32) 

The first term on the right-hand side is the subtraction term for final-state singularities 
with final-state spectators and has the same form as Eq. ( |6.10| ): 



/ dfff = '^il H H H f d(f)m{pi,---,Pij,---,Pk,---,Pm+l]Pa + Q) 



1 



XtT- ^!r\pi^ ■ ■ ■ ^Pij^ . . . . . . ,Pm+i;Pa) 

J{rn+1} 

X J [(iPiiPij,Pk)] 1^ij,kiPl, ■ ■ ■,Pm+i;Pa) , (6.33) 

where the flux factor fulfils the following scaling property: 

HVPa) = VHPa) ■ (6.34) 
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The integration of this term follows the steps of the previous section and leads to the 
analogous result, which, for later purposes, we write in a different form: 



m+l 



^^a = E / / \d(y^,{xpa) ® U""' 5(1 - X) 



(6.35) 



with /m(e) given in (|6] 



The second term on the right-hand side of Eq. (|6.32|) is the subtraction term for final- 
state singularities with initial-state spectators. We have seen in Eq. ( p. 46 ) that the phase- 
space factorization leads to a convolution over the boost parameter x, thus 



Mir, 



m+l 



i,3 



X ■ 



1) 



) Vij ) 



X / [dpi(pii;j5a,a;)]P,"j(pi, . . .,pm+\]Pa) ■ (6.36) 



Using the definition of the dipole function in Eq. ( ^.40 ), following the same manipulations 
as for the integral of the dcx^ term and finally making use of the scaling property in 
Eq. ( |6.34|) , we find that Eq. ( |6.36|) can be recast in the following form: 



dc^r = E / dx / \d(T^,{xpa) ® f^""' Im,a'{x- t 



/m+l - —JO 

The new insertion operator is 

as (47r) 



(6.37) 



Im,a'ix;e,iJ, ■,{pi,mi},pa) = - 



2n T(l 



2 \ <s 



^ Vj{x; Sja, rrij, {mp}; e) , 



(6.38) 

where sja = ^pjpa . We emphasize that pa is the original initial-state momentum of the 
incoming parton, while the final-state momentum pj belongs to the boosted phase space 
with incoming momentum xpa- The flavour kernel for a massive quark (or antiquark) j is 



Vg{x]Sja,mg,{};e) = Cp 
+5(1 - x) Cf 



JgQ[x,mq/y^s~j +6{l-x) Tq{y^,mq;e) + K, 



e Sja + 



(6.39) 



The function Vq^^Js^^mq] e) is that defined in Eq. ( |6.29| ) with the /i ^/sj^ substitution. 
Here and in the following, {} denotes the empty set. For a massless quark j = q the flavour 
kernel is given by 



Vq{x; s,a, 0, {}; e) = Cf [J^q{x, 0)J ^ + 6{1-x 

If j is a gluon, then 
Vg{x; Sja, 0, {mp}] e) = Ca 



CFfi-^l+r,(6) + ir, 



(6.40) 
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1 _ !L 



F=l 



niF 



+ Tg{{mF}]e) + Kg 
+ 5(x 



a;NS / ""^F 



QQ 



2 JV^" 2 



10 



(6.41) 



where x+ = 1 — 4m|,/sja. Note that the last hne, containing the difference of two 5 
functions, gives a contribution to Eq. ( |6.37| ) that can easily be calculated and has a smooth 
massless limit, since A^^" can be replaced by Np as mi? 0. 

The third term on the right-hand side of Eq. (|6.32|) is the subtraction term for initial- 
state singularities with final-state spectators. The relevant phase-space factorization, given 
in Eq. ( |5.46| ), is again a convolution. 



da 



m+l 



('^Zt \ E EE / dx /d0„(pi,...,Pj,...,p^+i;xpa + Q) 

^s(«)^(Pa) {^1} -^0 J 



X ■ 



S{m+1} 

X j[dp,{pj-pa,x)]Vf{pi,...,Pm+l]Pa)- (6.42) 

The evaluation of the integral over [dpj(pj; p^, a;)] follows that in the massless case, using 
the explicit results in Sect. 5.3.3| . Then the counting of the symmetry factors is also 
straightforward and using the scaling property in Eq. ( |6.34| ), we find that Eq. ( |6.42| ) can 
be recast in the following form: 



m+l 



dcr^'" = E / da; / \<la^,{xpa) ® Im,aa'{x]e) 



(6.43) 



where the insertion operator is 

Im,aa'{x-e,H^-{pi,mi},Pa) = Jl^^^ . E ' 

27r 1 (1 - ej ^ 

(6.44) 

The functions V"'" (x; Sja, rrif e) are related to the /"''(x; /i^; e) functions given in Eq. ( |5.88| ), 



s 



-rV'" (x; Sja,m,-e) . 
J- 



V'" (x; Sja,mj-e 



jaa 

3 



X 



rrij _ 



"JO. 



(6.45) 



Note that da^'" is the only contribution to da^ in ( |6.32| ) where the flavours of the partons 
a and a' can be different. However, both a and a' are massless, since we do not allow 
for massive partons in the initial state. The fact that V^'° does not receive contributions 
from the splitting g —>■ QQ into massive quarks Q is not a restriction of the presented 
formalism, it simply means that the logarithmic corrections ag In ttiq of this origin are not 
extracted from the real correction da^ in ( |2.2| ). These terms could be extracted by a further 



subtraction into suitably-defined parton distributions, as discussed at the end of Sect. 2.2 
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Next we sum Eqs. ( p.35|) , (|6.37| ), (|6.43|) and the well-known expression for the collinear- 
subtraction counterterm, 



, r/ 2^ "s (47r)' ^ , 



^al,{xva), (6.46) 



to obtain the sum of the full subtraction term and the collinear counterterm in Eq. (|6.31|) . 
In Eq. ( |6.46| ) [ip'vs, the factorization scale and the functions i^'^^ [x) define the factorization 



scheme. The MS scheme is defined by K'^i^ix) = 0. More on Eq. (|6.46|) can be found e.g. in 
Sect. 6 of Ref. |T^. We see that /^_,_]^dcT^(pa)+ Jm^cTaiPa) is obtained from the leading-order 
expression Jm^fJai^Pa) by replacing the leading-order matrix element squared 



ns{a) 



{...■,Pa\\...]Pa) m,a 



by 



^ ^ / ,\ m,a' {■ ■ ■ J -^Pal I (-^j ^) \- • ■ i •^Pa) m,a' 



7 ns[a' 



(6.47) 



(6.48) 



and performing the x-integration. Here the insertion operator [x; e) also depends on 
the renormalization and factorization scales, colour charges, momenta, masses and flavours 
of the QCD partons, 

i""'"'(x;e,/i^/i|; {pi,mi},pa) = 6""' 6{1 - x) i"m(e,/i^; {purrii}) 



(x;e,/i ■,{pi,mi},pa) + 1 m.aa' 

{x;e,fi ; {pi,mi},pa) 



«s (47r)^ 



-(^) P"'^'(x) + Krs'(:^) 



+ 0(e). 



(6.49) 



2tt r(l -e) 

One can understand the structure of the insertion operator, if we rewrite it as 

/"•"'(x;e,/x^/i|; {pi,mj},pa) = (5'"''(5(1 - x) /„+a(e, /i^; fe, mj},^^) (6.50) 
+ P''X'{x;fil-{p^},xpa)+K';^^'{x;{p,,mi},Pa) + 0{e) , 
where the operator Im+a{^, fJ''^', {Pi,^i},Pa) is similar to the operator defined in Eq. (|6.16|) , 
i"m+a(e,Ai^; {Pi,mi},pa) = Imie,n^; {Pi,mi}) + i"a(e,Ai^; {Pi,mi},Pa) , (6.51) 

where 

Ia{e,n^;{Pi,mi},pa) = fj^^T^ 



— T- T 

J 



( ~ ) ("'^i(-5ja,mj,0,{mi7};e, k) - ^ ) +T j{^,mj, {mF}]e) 



^30,, 



a 



+7j In — + 7j + Kj 

2 \ ^ / 2 \ 

V.(s,„0,m,-,{};e,2/3)-^ +^ 



Tl 



/i 

+7a In h 7a + -f^a 

Saj 



(6.52) 
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The only difference between and Im+a is that the latter depends on the additional 
initial-state parton a such that in the functions Va and the set {mp} is empty and 
K = 2/3. The singular terms V^^^ do not depend on the value of {mp} and k, and crossing 
the momentum of a parton from the final to the initial state does not change the singular 
terms in the one-loop QCD amplitudes, therefore, this insertion operator cancels all the 
singularities in the virtual contribution /^d(jj(pa). As a result, the operators P^" and 
are finite in four dimensions, as can also be seen from their explicit forms below. 



The operator, which is defined by 

a. 



27V 



P''^ i^)7f^T.TrTa' In 



XS 



(6.53) 



JO. 



is independent of the presence of massive partons and contains the full dependence on the 
factorization scale fip- In order to obtain this form of the fip dependence, we have used 
colour conservation (J2jTj = —Ta') and expanded in e. 



e 



XI 



P^''\^)^Y.TrTa' In ^ + 0(e). 

(6.54) 



The K'l^ operator contains all remaining terms, including the factorization-scheme 
dependence and further finite-mass corrections: 



Kr'(x;fe,m.},p,) = ^ ir'' 



\X 



K%,[x)- 



PZ (^) In ■ 



X^Sjd 



1 - x)sja + m] 



Sja - 2m j J Sja + m?,+ 2m] 
X I In ^— + 

Sja 



Tj ■ Ta' IC^'""' (x; Sja, ruj, {mp}) 



+ 7a 5'^'^' 5(1 -X) 

2m i 



Sja + mj + mj ^ 



. (6.55) 



Comparing the expressions for the various insertion operators on the right-hand side of 
Eq.( |6.49 ), we observe that the finiteness of K'^ is highly non-trivial: terms from the 
integration of all possible dipole contributions have to be combined, together with the 
collinear subtraction, before taking the e ^ limit. 

In order to simplify the final formulae, in Eq. (|6.55|) we used the functions P^^^{x) 
and P"-"-' {x), defined in Eqs. ( |5.89|) and ( |5.94|) , respectively, and also the K"'"' (x) functions 



defined in Ref. \\1 



K'"' ix) 



P;4'(x) ln^^+P'""'i 



X 



+5° 



X 

2 inl 



X 



1 — X 



X 



6{l-x)[^a + Ka--n^Tl 



(6.56) 



Recall that the functions P''*^(x) come from the e-dependence of the splitting functions 
and are given in Eq. ( 5.93| ). 
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The functions /Cg'"'(a;; sja, rrij, {mp}) depend on the flavours of partons a, a' and j. For 
the sake of clarity, we consider the various cases separately. If j is a quark, then /C^'"' does 
not depend on {mp} and we have 



/C^'''(x; Sja,mj) = 
Sja,mj) = 2 



(6.57) 



'ln(l — x) 



1 — X 



ln(2 - x) 

1 — X 



+ 



r X ^ 



+2 



\nJ^^-^^^-;^5{i-x) 



l-xj+ (2-x)sja + m^ Cf 



^^In 



Ca (1 - + mj ' 



A^g'^(x, Sjfl,, TTijr) ICg''^ (^X, Sja, TTlj) -\- JCg'^ (^X, Sja, TTlj^ 

Of 



(6.58) 
(6.59) 
(6.60) 



If j is a gluon, then 



X:«''^'(x;s,„0,{m^}) = -5'^'^'^ 



A F=l 



X 



m, 



- In 

3 I s 



+ 3, 



j: 



a:NS / "^-F 



+ 



5(1 -x)- 1 ^ 



(6.61) 



where x+ = 1 — 4m|,/sja and iV|f is the number of flavours with Sja > 4m|,. Note that in 
this case mj =0, so the terms proportional to 1/T^/ in Eq. ( |6.55|) vanish. In the zero-mass 
limit the functions JCj'"" simplify to 



/C"'" (x; Sja,0,{}) 



1 — X 



+ S{l~x) 



(6.62) 



Using Eq. ( |6.62| ), we see immediately that in the massless limit the operator -K'^" tends 
to the corresponding massless operator smoothly. 



x; {pi,mi],pa) 



(6.63) 



where K'^ (x) is the same operator as defined in Eq. (8.38) of Ref. [TT 



Finally, we explain the actual evaluation of the '+ '-distribution that contains the kine- 
matic variable Sja in more detail. As mentioned at the end of Sect. |5.2.3| , Sja has to be 
kept fixed during the x integration if Sja appears inside the '+ '-prescription, as is the case 
in JCj' for a = b. Schematically we have to evaluate an integral of the form 



^^dx J d$(x) [j[x,s%^)]^ \M{^{ 



(6.64) 
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where the label (x) of SjJ indicates that in caclulating this variable the final-state mo- 
mentum pj belongs to the phase space of the x-boosted frame, but pa is the original 
initial-state momentum of the incoming parton. At first sight this seems to be very incon- 
venient, but in practice the procedure can be simplified according to 

dx y d$(x) ds,aS{s,a-sfJ) [J{x,Sja)U 

= dx ds.a J {X, -S,a) [j d$(x) \M{^{X))\^ 6 [S,a - sfj) 

- /d$(l)|A^($(l))r5(s,,-.«)} 

= J\x [JdHx)j{x,sfJ)\Mmx))\'- ldmj{x,s^l^)\Mm))\']. (e.es) 

Note that the variable sja that appears in the function J{x, sja) is the one in the respective 
phase-space integral. A non-trivial change in the endpoint contribution occurs if this vari- 
able Sja is ehminated in favour of an x-dependent variable /(x), such as /(x) = — x Sja 

as described in App. ^ More precisely, if sf] in Eq. (|6.65|) is replaced by /(I), the endpoint 

changes; if sf] is replaced by /(x) the endpoint remains unchanged. 



6.4 Jet cross sections with two initial-state hadrons 



The derivation of the insertion operators Im+a+b, Pm+b and K^+b^ relevant for calculating 
the NLO correction of Eqs. ( |6.(j| ) and ( |6.9| ), does not contain any new features as compared 



to the combination of arguments given in Sect. 10 of Ref. |Tl| and in the previous subsection. 
Therefore, we simply list the corresponding results. 



The Im+a+b insertion operator is again very similar to Im of Eq. (|6.16 ). The only 
difference is that it depends on the additional initial-state partons a and h such that in the 
functions Va{b) and Ta(b) the set {mp} is empty and n = 2/3, which does not influence the 
singularity structure: 

Im+a+b{e,fi^;{pi,mi},pa,Pb) = 

Im{e,lJ^'^; {Pi,mi}) + /a(e,/x^; {pi,mi},Pa) + Ib{e,IJ''^; {Pi,mi},pb) 



as 



(4vr) 



2n T(l 



Sab. 



+ ia^ b)j . 
(6.66) 

The operator Pm+b is completely analogous to Pm of Eq. ( |6.53|) apart from the trivial 
dependence on the additional initial-state parton. 



Pmlbi^'^ f^F'^ {Pi},XPa,Pb) 



Pm (■^5 



fil;{p^},xpa) + ^P'^'{x)^Tb-Ta' In-^. (6.67) 



XSab 



Finally, the K^+b operator is 

^rntbi^'^ {Pi,mi},Pa,Pb) 



KZ\x\ {P^,mi},Pa) 



• T,j^P-'(x) ln(l - x) + 5'^ 



'ln(l 



1 — X 



. (6.68) 
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We see that K^+b has the form of of Eq. ( |6.55 ) with the additional term in the second 



hne containing parton-parton correlations between the initial-state partons. 



7 Summary 

In this paper, we have presented an extension of the dipole subtraction method for calculat- 
ing arbitrary (phenomenologically relevant) jet cross sections at NLO accuracy in arbitrary 
scattering processes involving heavy partons in the final state. In the case of lepton colli- 
sions, we have set up the formalism such that for those cross sections for which the mass 
of the heavy parton does not set the hard-scattering scale, i.e. the massless limit is IR safe, 
we can simply set the masses of the partons to zero to recover the massless limit of the 



computations, discussed in great detail in Ref. ||TT[. Thus the implementation of our for- 
malism into a general purpose NLO partonic Monte Carlo program leads to a code that is 
smooth in the massless limit or, more importantly, numerically stable for any values of the 
hard-scattering scale. In hadron collisions, similar behaviour can be achieved by matching 
the partonic calculation with a suitable definition of the heavy-parton distributions in the 
massless limit. This feature will be particularly important at the future colliders (LHC, 
NLC), where the ratios of parton masses to other relevant kinematical invariants can run 
over a very wide range of values. 

The factorization of the QCD matrix elements on soft poles in the presence of massive 
partons has been known and used for a long time. In those kinematical situations in 
which two partons become collinear and at least one of them is massive, the collinear 
divergences in the matrix elements are screened by the finite value of the parton mass and 
the regularization of the real corrections in these regions is not, strictly speaking, necessary. 
Nevertheless, the cross section receives a logarithmically enhanced contribution from these 
phase space regions when the parton masses become small compared to other relevant 
kinematical invariants. In order to assure the smooth massless behaviour of the cross 
section, we introduced the notion of quasi-coUinear limit and presented the factorization 
of the QCD matrix elements on quasi-coUinear poles (factors that become real poles in 
the massless limit). The smooth massless limit also serves as a powerful check of both the 
analytic calculations and any numerical implementation. 



Our extension is formulated very closely along the lines of Ref. [rT|. We changed the 
formalism to the least extent that was necessary to incorporate the mass corrections. Tech- 
nically the generalization is cumbersome, but the complications only concern the derivation 
of the insertion operators (done in this paper) and the actual application is not much more 
involved than it is in the massless case. As a result, an existing general purpose NLO 
Monte Carlo program, such as NLOJET++[T^, IS], written for computations in the massless 
theory, can be changed straightforwardly to incorporate parton masses. In particular, the 
general form of the dipole subtraction terms remains unchanged and the dipole splitting 
functions receive trivial mass corrections. Only the definition of the emitter and spectator 
momenta, needed for writing the matrix elements and jet functions of the subtraction term 
change in a somewhat cumbersome, but straightforward, way. 

The NLO contribution containing the virtual corrections and the insertion operators also 
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has the same form as in the massless case. We have presented these operators exphcitly, 
thus, for constructing a numerical program to calculate the NLO corrections to arbitrary 
jet quantities in a given process, the only additional ingredients that we need are simply 



related to the evaluation of the original matrix elements, as listed in Sect. p73 
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Appendix 



A The eikonal integral 

In Sect. ^.1.3| we have defined the eikonal integral P^^{nj,Hk) in Eq. ( |5.26| ), but did not 
give the general expression, since only the part symmetric in /ij ^ fik is needed in the final 
result. An intermediate result of the general integral, which is given in the following, is 



2e 



TT 

12' 



+ Lia (1 - (fij + ^ikf) - Li2 ( 1 - pf) - 2/(/ij, pk 



0( 



(A.l) 



with pj, pfc, p defined in Eq. (|5.3CI|) and f{pj,pk) is an auxiliary integral 



t_ = (/i, + pkY . (A.2) 



t + p]-pl + ^\{t,pj,pl) 
The explicit result for f{pj,pk) reads 



l2 



Pjp) 

(A.3) 



where the massless limit pj — > yields f{0,pk) 



TT 



76. 



The symmetric part of the integral defined in Eq. ( [A.2|) is much simpler than the full 
result and is given by 



^ [/(Pi, Pfc) + /(Pfc, Pj)] = ^ Li2 (l - {pj + pkf) - Li2(-p) + Li2(l - p) - ^ , (A.4) 
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with p as defined in Eq. (|5.30| ). We used the fact that Vij^k 
is a massless parton, as in the case of the eikonal integraL 



Vik,j if the emitted parton 



B Alternative integrated dipole functions 



At the end of Sect. |5.2.3| we have emphasized that the forms of the J^^ endpoint parts of 
the x-distributions depend on the convention for which kinematical invariant is kept 
fixed during the x-integration which is performed over a plus distribution [/(x, -P^)] + . This 
becomes obvious if we substitute the invariant = 2paPij used in Sect. |5.2.3| with another 
invariant such as, for instance, 



p2 = g2 ^ (Pij - xpaf = - 2xpaPij , 



(B.i: 



which was used in Ref. [^. In this appendix we provide the necessary formulae for this 
particular variant of treating the integrated dipole functions. 

The transition to the parametrization in terms of the new = requires some 



changes in the separation of the one-particle phase space described in Sect. |5.2.1| . The 
rescaled parton masses have to be normalized to P^, so we define 



The previously defined rescaled masses /i„ are, thus, replaced by yU„ using 



2 



— 2 
1 + n% 



(B.2) 



(B.3) 



in Eqs. ( |5.48| ) and ( p.49| ), leading to 



X 







— e 




/ d'^-^n / d2i 


Z+{x) - Zi 




Zi - z_(x) 


1 Jz^ix) 







From these expressions the upper limit x+ for Xij^a is derived in terms of fin, 



1 + (/ij + lljf 



(B.5) 



The dipole functions V^, defined in Sect. ^.2.2| , remain unchanged. The integrated 
dipole functions 



[dp,{Q]Pa,x)] 



1 



rat 



^^^^■^ - 2I r(i - e) (m?. - g2 ) ^^^(^^ 



(B.6) 
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are given in terms of new auxiliary functions Iij{x; e) and their calculation proceeds along 
the same lines as described in Sect. p.2.3| . The endpoints are split off according to 



Pqq{x- e) = Tr ([J^q(x, /xq)],^ + 5{x^ - x) [j^liM. e) + JqI^^q)]) + 0(e) , (B.8) 
where a;+ = 1/(1 + 4/iQ). The continuum part for the Q gQ splitting, 



1 — X 



In 



+ 



^^/2 + 2fil-x 



+ 



:i + -^ini-x) 

2{l + fil-xy 



l-xj+ \ l + l^tj J 
and the singular part of the corresponding endpoint contribution, 



1 1 



TT 



ilnfi±^U± 



2e 6 



:7g + 



(B.9) 



(B.IO) 



result from their counterparts [Jggix, ^q)]+ and J^g {fig; e), given in Eqs. (|5.58|) and (|5.59|) , 
upon the simple substitutions fig xJiq/{1 + Jiq) and /Iq — >• fig/ {I + ftq), respectively. 
The non-singular endpoint contribution is given by 



JgO (/^q) = — - 2 Ll2 



1 



-2Li2 



1 + 2/i! 



--In^ 



1 + 2/ii 



+ 1.^(2.,^) in 

a;NS 



:B.ir 



Note that this does not follow from J^q (/ig), given in Eq. ( |5.60| ), upon substituting /iq, 
but a finite difference remains 



'gQ 



/ 1 
2 Lia . 

3 M 1 , r.2 



1 + /^Q. 



1 



+ o/^Q(2 + 4)ln 



1 + 

1+4 



2Li2 



lln 

2 



,1 + 24. 



r,2 



1 + 2/x^ 



(B.12) 



For the g — * QQ splitting the same features are observed. The continuum part, 

= ^ I — ^wi 



Axfi'^ 



(l-x)2 
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Q 

1 — X 



[B.13) 



x+ 



results from [Jqq(x, AiQ)]x+, given in Eq. (|5.62|) , upon substituting fig x/ig and 
The singular endpoint part, 



3e 



10 
"9"' 



(B.14) 



is obtained from JggifJ'Q', e), given in Eq. ( p.63|) , by the replacement /ig — > x+fig. However, 



the non-singular endpoint part differs from J^Q^^fiq), given in Eq. ( ^.64| ), in a non-trivial 
way: 



4 



,ln 



1 + 4/i^ + l 



. (B.15) 



The difference between JqQ^ifJ'Q) and the reparametrized JgQ^ifJ^q) reads 



ta;NS / - 



'QQ 



QQ 



8 



10 



+ 



Ta;NS 
QQ 

10 



f^Q 



1 + 4/i^ 



4 

+ 3 



2/i 



r,2 ^ 



16 _2 



1 + 44-1 



In 



1 + 4^1 

2flQ 



:B.16) 



Since the non-trivial differences between the two endpoint parametrizations appear only 
in the integrated dipole functions for final-state emitter and initial-state spectator, the 
necessary changes in Sects. and |6.4| can be easily inferred. The new functions 
follow from their previously defined counterparts /CJ* by simple substitutions and adding 
the extra terms AJ^q^ and AJ"''^^ 



'QQ 



'ja—im^j-Q%)/x 
JS ( 

gQ 



rrij 



ICf{x-Ql,0,{mF}) = /Cf(x;v,0,{m^}) 

rji Np 



^ /x, 

la' ' " 



rap 



F=l 



'Q ja I 



:B.17) 



:B.181 



where = ""^j ~ xs^a-, not to be confused with = Sja + ^| used in Sect. of the 
main text. Note that in the scheme considered in this appendix, all masses of quark are 
accessible at all values of so Np appears in Eq. ( [B.18| ) rather than N'p . In fact all 
occurrences of Np!" in / and K should be replaced by Np in this scheme. 
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C Dipole functions in SUSY QCD 



In this appendix we define the dipole functions Vij^k and V^j needed to construct the dipole 
subtraction function specified in Eq. (|5.1|) for SUSY QCD calculations (see for example 
Refs. 



m 
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The new dipoles fall into two classes: one that involves gluinos and 
gluons and the other that involves squarks and gluons. We consider only those dipoles that 
lead to 1/e poles after integration in d = 4 — 2e dimensions. Owing to the large masses 
of the SUSY particles in typical supersymmetric extensions of the SM, the massless limit 
is not too relevant even at LHC or NLC energies. Thus we do not define dipole functions 
corresponding to either g gg, qq or splittings involving a squark, a gluino and a quark. 

In principle, to calculate one-loop amplitudes in supersymmetric theories one has to 
use a regularization prescription that respects the supersymmetric Ward identities, for in- 
stance, dimensional reduction (DR). Within the dipole formalism, the use of different regu- 
larization prescriptions only affects the actual calculation of the cross section contribution 



^NLO{m} ( |6.8| ). In this paper we derived our formulae using conventional dimensional 

regularization (CDR), therefore, the one-loop amplitudes \J^m,ab{{Pi,^i}]Pa,Pb)\fi-.ioop) 
are also needed in CDR. To derive the insertion operator /(e) for SUSY QCD one has 
two options, either to use DR throughout and the transition rules between CDR and 
DR elucidated in Refs. |Q [Tl|, 0, or to use CDR throughout and correct for the 
violation of the SUSY Ward identities in the calculation of the one-loop contribution 
\Mm,ab{{Pi^^i}iPa,Pb)\fi-ioop) explained in Refs. |5|, || 



We choose to do the latter. 



As discussed in Sect. H, we do not consider massive partons in the initial state, therefore, 
new dipoles of the P^* type do not appear in the supersymmetric theory. The dipole 
functions that involve gluinos are obtained immediately by changing the colour factor Cp 
in Eqs. (|5.16|) and ( |5.50|) to Ca- As a result the integrated dipole functions are also obtained 
by this simple change in Eqs. ( 5.23| ) and ( |5.56 ). 



The flavour and the spin of the spectators do not influence the dipole functions, there- 
fore, we have to consider only flnal-state emitter squarks as new cases. Thus, we deflne 
the splitting functions Vjj ^ and Vf^- for the q ^ gq splitting. The case q ^ gq is formally 
identical to g ^ gq. 

When the spectator is also in the flnal state, then we deflne 



{sly gg^k\s') = Svr/i 'OsCf i r - 

= ^ gq,k)^ss> ■ 




(C.l) 



We write the integral of {Vgq,k) / i'^PgPg) over the one-parton subspace in Eq. (|5.11| ) in the 
form of Eqs. ( |5.22| ) and ( p.23| ). The eikonal integral depends only on the mass of the 
emitter, but not on its flavour or spin. The coUinear integral is obtained as 



-2e 



2/ir^^ + 6-21n {1 - - 



+ P^ + Oie). (C.2) 
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If the spectator is in the initial state, we define 



2 ^ij,o> 



PiPj J 



(C.3) 



We write the integral of (V^^) / {2pgPg) over the one-parton subspace in Eq. ( |5.48| ) in the 
form of Eqs. (|5.53| ) and (|5.56|) . The continuum part is 



1 — X 



1 + ln(l - X + /i? 



+ 



1 — X 



ln(2 + /x^- - x) . (C.4) 



The endpoint parts are 



where 



Ta;NS/ N ^ oT- 

Jgq (Z^?) = y - 2 Li2 



7g = 2Cf 



2Li2(-/i?-)--ln2(l + /i?), 



TT 



F • 



(C.6) 
(C.7) 



Using these results, we obtain the functions needed for the construction of the insertion 
operators and K""'^ . The singular functions V^^-* are independent of the fiavour and 
spin. The non-singular function V,-'^^'' for the gluino is 



Vf^)(.,,,m„m,) = Vr)(.,.,m„m,) 



(C.8) 



and for the squark it is 



Vf^^ [Sjk, rrij, mfc) = In ^ - 2 In 

q ^jk 



+ 



jk 



Sjk 



1 

+ 

Vjk 
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\np' ln(l + p') + 2 U,{p') - Li2(l - pj) - Li^ll - P,) - 



(C.9) 



if k is massive, and 



2 In 



Sjk 



jk, 



jk 



(C.IO) 



if k is massless. The singular function Tj for the gluino is 



(1 m 
e p 



2 +7,ln-f = C, 



p- 



1 I , 

- + -ln^-2 

e 2 p^ 



(C.ll) 



and the fiavour constants are 



75 = 



(C.12) 



53 



The singular function Tj for the squark is 
rg(/x, Trig; e) = T? 



mj: \ m 



- 2 



Finally, the functions /C"'* for the gluino are 



and those for the squark are 



2[sja{l-x) +mj] 



212 



(C.14) 



(C.15) 



D Example Applications 

We illustrate our method with the three simplest examples, e~^e~ annihilation to heavy 
quarks, e+e" — > QQ, to heavy quarks and a jet, e~^e~ QQg and heavy quark production 
in deep inelastic scattering, ep — ^ eQQ + X. 



D.l e+e- QQ 

We begin by recalling the tree-level result to set the notation. Since our aim is to illustrate 
our method, we try to keep the analytical formulae simple by taking the average over 
event plane orientations and neglecting electron polarization. The extension to oriented 
and polarized observables is straightforward. 

In order to match consistently with the one-loop corrections that are available in the 
literature, it is necessary to evaluate the tree-level matrix element in d dimensions. In fact, 
since all singularities can be cancelled before averaging over event orientation, it is sufficient 
to consistently evaluate the hadronic tensor in d dimensions. We are then free to take the 
rest of the process, as well as the remaining angular integrations, in 4 dimensions. Since 
after averaging over event orientation and polarization we only encounter Dirac traces with 
no or two 7^ matrices, we can easily eliminate 7^ by taking it to be totally anticommuting 
in d dimensions and using 7I = 1. 

With this prescription, the (i- dimensional matrix element is identical to the 4-dimen- 
sional one which, labelling the momenta by e+ -|- e" — > 7*/Z(g) — >• Q{pi) + Q{p2), is given 
by 15|] 

\M,\' = (^?^^ + (1 + 2^1) - (6/.y , (D.l) 
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where /xg = itiq/^^ with s = g^, and the normahzation is such that the cross section is 
given by 



0" 



LO 



(Tov\M2\^fP{p,,P2 



(D.2) 



where v = y 1 — is the velocity of the heavy quark in the centre-of-mass frame and 
Fj is our (infrared-safe) observable. The coupling constants and point-like cross section 
appearing in Eqs. ( |D.1D and ( |D.2| ) are given by 



= Ql-2g:g?QQRe{x{s)} + i9f + 9f)gf\xis)\^ 
9^^ = {9f + 9f)9f\x{s)\^ 



3s 



:d.3) 
:d.4) 

:d.5) 



where g{, = T/ — 2Qf sin^ 6yu, g[ = , Qf is the electric charge of fermion type / and 6^, 
is the weak mixing angle. The function parametrizes the Z° propagator and coupling 
factors, 

1 s 

^(^) = A ■ 2n ^ T^' 

4 sm tly^ cos^ tiyj s — rriz + inizi- z 
with mz and Tz the mass and width of the Z°. 

The NLO real emission process is e+e" — »• 7*/Z(g) — > Q{pi) + Q{p2) + 9{P3), with 
matrix element M3{pi,P2,P3)- In discussing the NLO corrections to the QQg process, we 
will need to discuss the c/- dimensional behaviour of A^3 in more detail, so we explicitly 
indicate the number of dimensions in which it is evaluated. In four dimensions we can 
rewrite the known result 



as 



\M^\pi,P2,P3)\'' = Cp' 



x\ \M2\ 



ma. 



2(1 - 



(1 — Xi) \2 — Xi — X2 



Xi 



+ (Xi X2) 



+ 



1 — X2 1 — Xi 
1 — Xi 1 — X2 



+ 9^^^f^Q 



(D.7) 



where Xi = 2pi-q/q^. In terms of these variables, the phase space is given by 



d$(3) 



-dxi (1x2 6(a;+ - X2)Q{x2 - x_) 0(1 - xi)Q{xi - 2^q), 



with 



(2 - xi)(l -xi + 2fil) ± (1 - xi)Jxl - Aiil 



2(l-xi+/i^) 
The three-parton cross section is then 

= (To / d$(=^) \M,\^ Ff\p,,p2,P3). 



(D.9) 



(D.IO) 



where we used a different normalization from that in Eq. (p 
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According to the dipole subtraction method, this integral is rendered finite by subtract- 
ing from it an auxihary cross section constructed from the two dipole contributions, ^31,2 
and ^32,1, defined by Eq. (|5.2|) . The colour and spin algebra are trivial in this case, and 
we obtain ^ 

2^31,2(P1,P2,P3) = (Vg3Qi,2) |-M2|^ 



with 



2p3-Pi 



^31,2 

1/31,2) ^^31,2 



1 + 5i + 



rrir 



P3-Pi 



Inserting the definitions of 2/31,2, -^i, "^31,2 and ^31,2 from Eqs. ( |5.12D , ( |5.14D and 
obtain 



(D.ll) 

(D.12) 

, we 



^^31,2(^1,^2,^3) = Cf' 



1 2(1 - 2/i2 



1 - X2 1 2 - xi - X2 \ 



1 - 4/i| X2 - 2/i| 



2 + — —-^ + 



X2 - 2/i5 I - X2 



(D.13) 



The associated dipole kinematics are given by 



^^=2^ + /.2_..,2 



xi - Afii 



P2 - 2^2?'' 



P31 



P2 - 2^2?^ 



The dipole 1^32,1 is obtained from 1^31,2 by the replacement pi ^ p2- 



(D.14) 



Combining the expressions for the real and auxiliary cross sections, we obtain the three- 
parton integral. 



NLO{3} 



(Jo Cf 
x||A<2|' 



1 



1 



dxi dx2 

/ 2(1 -2^1) 

1 — Xi ^2 — Xi — X2 

2(1 - 2fil) 



Ff\pi,P2,Pz) 



1 — Xi 



1 - xi - 2/i^ / ^ X2 



2 + 



1 — xi I 2 — xi — X2 \ a^i — 4/iQ 1 — 2/iQ \y ' xi — 2/i^ 

xFj^^(pi,P32) + (Xi ^ X2,Pl ^P31,P32 P2) 



1 — X2 1 — Xi 
1 — Xi 1 — X2 



Ff\p,,p2,P3)\.iB.15) 



1 — Xi 



It is straightforward to check that the soft singularity cancels between the different terms in 
Eq. ( p.l5| ), provided the observable is infrared safe (implying that F 



(3) 



Fj^^ in the soft 



(3) 



.(2) 



limit). Furthermore, for any quasi-collinear-safe observable (implying also that F 
in the quasi-collinear limit), Eq. (p.l5|) is finite in the small-mass limit. In fact, taking the 
quark mass to zero, one exactly recovers the massless result in Eq. (D.7) of Ref. JTT|. 

Next we have to evaluate the insertion operator /(e), which gives the integral of the 
auxiliary cross section and combine it with the virtual cross section. The one-loop matrix 
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element was calculated in Ref. |54] and is given byQ, 



|-^2|^l-loop) 



rUn 



Yd 



2Re{/i} \M2\ 



The form factors appearing in Eq. ( p.l6[ ) are given by 



(D.16) 



Re{/i} 



a 



1 + 



l + v\ 1 
m — 



2v 



Li2 



Re{/2} 



a^l — V 



In 



1-v 
1 + v 
1-v 



+ 



In' 



1-v 

1 + v 



2v 



+ In 



1 + v 



1-v 



In 



2v 



1+v 1+v 



2tt 2v 1 + v 
For later convenience, we define a function ff, 



(D.17) 
(D.18) 



1 + 
2t; 



In 



1 + w 



3 1-v 27c^ _ , 

2 - - In + In 

2 1+f 3 2 1 



1, ,1-v 



such that ff{v) vanishes in the small-mass (i.e. f — ^ 1) limit. 



(D.19) 



The general expression for /(e) is given in Eq. ( 6.16 ). In our case, the colour structure 
is trivial and we obtain 

«s (47r)^ 



2(l,2|J2(e,/^';te,m,})|l,2)2 = {M^l' x 2Cf 



271 r(l -e) 



TT 



— ) ( Vg(si2, mg, mg; e) - ^ ) + 7^rg(/i, mg; e) + ^ In + 5 - ^ 



TT 



Sl2 



Sl2 



(D.20) 



where Si2 = 2pi-p2 = s — 2mQ and we have suppressed the dependence on {mp} and k, 
which do not enter our cross section. The function Vq is decomposed into singular and 
non-singular parts according to Eq. ( |6.18| ). The e-expansion of the singular term is given 
in Eq. (|6.20| ). In our case it is 



1 + v' 
2v 



V 



e 1 + v 



-In^ 



1-v 

1 + v 
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In In 



1 + v 1 + v' 



The non-singular term is given in Eq. ( |6.21| ), and is given by 



(D.21) 



1 + v' 
2v 



21n- 



V 2(1+V) 
In — + 2 Li2 



1 + v il + v) 



+ ln(l - Wl-v^) - 2 ln(l - VT 



1-v 
1 + v 
1-v^ 
1 + v^ 



2Li2 



2v 
1 + v 



~6 



In 



2-VT^ 



+ 2 



1 _ ^2 _ 



+ 



TT 



(D.22) 



1 + v'^ 

**Note that we have inserted a factor of (47r/x^/TOg)^/r(l — e), where ^ is the dimensional-regularization 
scale, which is necessary for consistency with our notation. 
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Note that in the massless hmit, V^'"^^-' vanishes. The function Tg appearing in Eq. ( D.2C ) 
is simply given by 



Tg{fi,mQ;e) = Cf 



i + ^ln!^ 
e 2 



(D.23) 



Combining the virtual and auxiliary cross sections, we obtain a two-parton cross section 
that is finite as e — 0. Setting e = 0, we obtain (recall that the leading order cross section, 
0"'"°, is defined by the observable Fj^\pi,p2)), 



a 



NLO{2} 



-Vf TT^ 3, 2(1 + V^) 

In TTT 

12 2 (1 + v)^ 



+ In ^— I ^ . ^ In ^— + In ) + ff{v) + V^^^' 



1 + v 



1 + v 



+2Re{h}vao({g'''' + g^^ 



1 + V^ 

5 



4/i^))^f (Pi,P2). (D.24) 



Note that not only have all singularities cancelled, but also all terms that are singular in 
the small- mass limit. Furthermore, in this small- mass limit a^^'^^'^^ agrees with the exactly 
massless prediction given in Eq. (D.ll) of Ref. |TT|, 



lima^LO{2}^^LO^ 



TT 



(D.25) 



As a check of our results, we can use them to obtain the total cross section, by setting 



.(3) 



1. The integral in Eq. ( |D.15| ) can be performed analytically, but the result 



does not have a compact form, so we give its expansion in powers of fig, 



a 



NLO{3} 



-.VV 



AA 



+ ^i'q [-2 In /i^ + 3) + 4/i$ + 0(/i^ In /i^) 



+9 



AA 



f^l (-81n/iJ-15) + 0(/i^Qln/.y 



(D.26) 



It is worth noting that no terms linear in fiq arise. Although there are no such terms in 
the real or virtual cross sections, they could in principle arise in the auxiliary cross section, 
cancelling between the two- and three-parton cross sections. If this did happen it would 
worsen the numerical convergence in the limit of small but finite masses. Such terms do 
appear for all higher odd- integer powers of fiQ, but are not problematic. 



Combining with a^^*^^"^^ from Eq. ( p.24|) , we obtain 



cr 



NLO 



- + 18^| + 0(/i^ln/i|) 



+9 



AA 



/i^ (-181n/.^-27)+0(/x^ln/i^) 



(D.27) 



in agreement with the result given in, for example, Ref. [Q. In fact, expanding the two 
results to arbitrary order in fig we obtain perfect agreement. 
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The example application discussed in this appendix has also been considered in Ref. 
within the context of the methods of Refs. [0 and As mentioned in Sect. the 

formalisms of Refs. 0, are not aimed at smoothly controlling the small- mass limit. 



This is explicitly shown, for instance, by comparing our Eq. ( p.24| ) with the analogous 
equations (the sum of Eqs. (15) and (16) for the phase-space slicing method of Ref. 0; 
the sum of Eqs. (16) and (24) for the dipole method of Ref. ||3^) in Ref. [56|, which are 
singular in the small- mass limit. 



D.2 e+e- QQg 

We continue to average over event orientation and polarization. As in the e~^e~ QQ 
case, we require the tree level matrix element evaluated in d dimensions, which is given by 

|Mf(Pi,P2,Ps)l' = \Mi'\puP2.P:,)?-i\M<i\p„p,,p,)\\ (D.28) 

At NLO, three different real-emission processes contribute: (a) e~^e~ —>■ Q{pi) + Q{p2) + 
9iP3) + giPi); (b) e+e" Q{pi) + Q{p2) + qips) + qipi); and (c) e+e" Q{pi) + Q{p2) + 
QiPs) + QiP^)- The quark q could be any flavour other than Q, massless or massive. Since 
the matrix elements, for these processes are rather lengthy, it is not feasible to explicitly 
show the cancellation of the soft and coUinear poles between the real-emission matrix 
elements and the auxiliary cross sections constructed from various dipole contributions. 
Therefore we only give the parts of the auxiliary cross sections, but do not reproduce A^4, 
which can be be found in Ref. IE? . 



For process (a), we have to evaluate ten dipole contributions, 1^31,2, ^^31,4, ^^41,2, ^^41,3, 
^^32,1, "^32,4, ^^42,1, "^42,3, ^^34,1 and ^34,2- The associated colour algebra is again straight- 
forward because the different colour projections of the three-parton matrix element fully 
factorize (see Appendix A of Ref. |11[). Thus we do not need to calculate any colour- 



correlated tree amphtudes and we obtain 

^^&(Pl,P2,P3,P4) = TT^fl-T^) (V33Qi,2) |A^3(P31,P2,P4)|', (D.30) 

1 Ca 

^^3il4(Pi,P2,P3,P4) = 7^ 7^ C^gsQiA) l^afei, P2, P4) (D.31) 

1 1 

^^34ll(PbP2,P3,P4) = - {n\Vg,g^^i\u) r^„{pi,p2,P3A). (D.32) 

^P3'P4 ^ 

The dipole contributions 1^32,1, "^32,4 and ^34,2 can be obtained from T'31,25 "^31,4 and 1^34,1 
respectively by the replacement pi P2- Likewise 1^41,2 and ^41,3 can be obtained from 
1^31,2 and 1^31,4 respectively by the replacement ps <-> p^. Finally, ©42,1 and 1^42,3 can be 
obtained from 1^31,2 and 1^31,4 respectively by the replacement pi p2 and ^3 ^ ^4. 



The splitting functions (V^gQ^ 2) and (V^gQ^ 4) are given by Eq. ( ^.16|) with = ttiq 



and rrik = respectively and (/^iV^g^^ i|z/) by Eq. ( p.l9| ) with rrik = mq. The tensor T^^y is 
given below. 
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For process (b) we have to calculate at most four dipole contributions, ^^34,1 and 1^34,2 
and, perhaps, ^12,3 and T'i2,4- We obtain 



^^34!l(Pl)P2,P3,P4) 



(D.33) 



The dipole contribution 1^34,2 can be obtained from 1^34 j by the replacement pi P2- 
The splitting function (/^iV^gg^ is given by Eq. (|5.17|) with = mg. Note that we 
have explicitly kept irtq non-zero in Eq. (p.33|) . Since all our results tend smoothly to the 
massless results in the small-mass limit, we can easily replace by zero if necessary. 

The kinematics for the dipole contributions ©12,3 and I^i2,4 replace the heavy quark and 
antiquark momenta pi and p2 by a massless gluon with momentum pi2. If our observable, 
Fj, requires the presence of heavy quarks in the final state then it will get zero contribution 
from these dipole contributions, leaving uncancelled logarithms of /Iq in the 4-parton inte- 
gral. Even in this case, it may be helpful to replace the observable by a pseudo-observable 
in which the subtraction qqg configurations do contribute, improving the numerical con- 
vergence of the three-parton integral and allowing the logarithms to be isolated in the 
pseudo-coUinear region. This may be more amenable to analytical treatment, allowing the 
logarithms to be summed to all orders ||5^. If the observable sums over all final states, then 



Pi2,3 and Pi2,4 are needed. If necessary, they can be obtained from T'34,1 and T'34,2 by the 
replacements Pi Ps, P2 Vi and mq ^ nig. 

For process (c) we have to calculate eight dipole contributions, 1^12,3, ^^12,4, ^^14,2, ^^14,3, 
^^23,1, 1^23,4, "^34,1 and P34,2- They are identical to those for process (b), with nig replaced 
by tuq. Note however that the external factors in Eq. ( |6.4| ) introduce an extra factor of | 
coming from the counting of two pairs of identical particles in the final state. The other 
dipoles can be obtained by appropriate permutations of momenta. 



The tensor T^i, appearing in Eqs. (|D.32| , [1X33D is the squared amplitude for the LO 
process e~^e~ — > QQg not summed over the gluon polarization. It is normalized so that 
—g^'^T^^ = IAI3P. Again averaging over event orientation and polarization and neglecting 
terms that are antisymmetric in n, u, which cannot contribute to the cross section, we 
obtain 

SttOs [9 +9 +9 ^^Q^^,u 

■ltiu{Pl,P2,P3j = -<-^F 



(l-a;i)(l-a;2) 



(D.34) 



with 



r 



VV 



AA-VV 

fJ,U 



2(l + 2/i^ 



P1P2 ^ P2K 



- 2- 



X 



1- X2 



-2 



1 - X2 

, 1 

1 — Xi 



1+2/i: 



I - xi - X2 + xl- 2/iQ(xi - X2) 
I - X2 

1 — X2 — Xi + xl ~ 2/iQ(x2 — Xi] 



1 — Xi 



Ml/ Ml/ 

P1P3 ^ pm_ 



P2P3 ^ P3P2 



+ (1 + Ixl + \x\ - xi - X2) g^"" + 4/i; 



2 P3P3 



-12 



P1P2 ^ P2P1 



+ 12!^^ + 12! 



X2 P2P2 



X2 



Xi 



2 ^P2P2 
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(D.36) 



Next we have to evaluate the insertion operator /(e), which gives the integral of the 
auxiliary cross section, and combine it with the virtual cross section, which was calculated 
in Refs. ^ . We use the notation of Ref . |^ . They obtain 



|A^3(Pl,P2,P3)la 



(1-loop) 



2N 



f 



2 1 

1 1 



17 



+ 2 In 



+ In 



(1 - Xi)(l - X2) 



-lE 



- — - 2/3 - (1 + In 



3iVe m (3 



a. 



1 



+ ^N, -\Mi;\pi,P2,P3)\^ + i^finite(Pl,P2,P3) 

Zn e 



\Mt\puP2,P3)\' 

(D.37) 



where Nf is the number of massless flavours and (3=^1 — 4/iQ/ (xi + X2 — 1) is the velocity 
of pi in the P1+P2 frame. The function -Fgnito contains all of the finite (as e — 0) terms. 



finite 



x2[g 



VV T7VV 



rpv V I AA rpAA \ 
-^finite "T y -^finite ) 5 



where 



:pVV/AA 
finite 



counter,finite,yV/ylA T-iext.,finite,y V/AA 



a. 



lc,VV/AA 



•sc,VV/AA 
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, (D.39) 



with i.-unter,finite,yyMA^ ^ext.,finite,yy/AA^ ^l^cVV/AA pScyV/AA ^^^^^^ EqS_(4.11), (A.2) 

and (A. 3) of Ref. In Eqs. ( p.37| - [n.39D the expansion parameter is a^^, i.e. charge 

renormalization is carried out in the MS subtraction scheme. 



In checking the quasi-collinear limit of our formulae, it is useful to note the logarithmic 
terms in Ffinite. We obtain^l 



Ffinite(Pl,P2,P3) = ^IM^^ \' l^lu ^ (^^C A - f + Np)^ 



--Ca 7£; - In 



2Cf — Ca , 
H m 



l-V, 



QQ 



lE - In 



'QQ 



+ 7£; - In 



47r/i2 



11 



sm, 



Q 



2Cf + — CA--TRiV/ + CAln. 
^ 6 ^Qg^Qg, 



Nf 



+Cf In^ - In /i J + -Tr ^ In - 

6 S 



a. 



+ ^\M^^\'\-CAhE-ln 



+ 2CFln/x^| + 



(D.40) 



ttNote that the result in Ref. [|9| is for a single heavy flavour. We have trivially modified it to incorporate 
Np heavy flavours, which simply involves changing the logarithmic term in from ln(zs//L(^) 

toE^:iln(m|/M'). 
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where the elhpsis represents terms that stay finite or vanish for /ig — and m?p/s 0. 
The first fine comes from the renormahzation of a^, with the renormahzation scale, and 
will remain in the final result. All the other occurrences of ^ come from the trivial fact 
that in Ref. the natural overall factor of (47r/i^/s)'^/r(l — e) is expanded in e and are 
cancelled by our insertion operator, I{e). All the remaining mass logarithms are related to 
the quasi-collinear limit and should also be cancelled by /(e). 

The general expression for /(e) is given in Eq. (|6.16| ). The colour structure again 
factorizes and we obtain 



3(1, 2, 3|/3(e, /i^ te, m,}) 1 1, 2, 8)3 = \Mf\^ x ^ ^^''^ 



2 \ f 



+Ca 



Itx rfl - e 
2 \ 

III. 

F — Ca 



"^Q,0) 



2Cf — C'a, .(NS) 



Iq In 



Iq In 



2 

+_^^^lni^ + i^^ln^+(l 
2Cf 2 

-(2Cf + Ca)— + 2rQ(/i,mQ; e) + 27, + 2ir, + rg({mi.}; e) + 7^ + 



Ca_ 
2Cf 

Ca_ 
2Cv 



a ( ' 
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where Sjk = 2pj-pk, V*^^^ is defined in Eq. (|6.20| ), Vj^^^ are defined in Eqs. ( |6.21| - [6l26| ), 
are defined in Eqs. ( |6.27| - |OUD and 7^ and Kj are defined in Eqs. ( |5.91D and ( p. 17] ). 

It is worth stressing that Eq. (p.41|) is the insertion operator for the QQg final state. 
It does not therefore account for terms coming from the splitting process g —>■ QQ in 
process (b) above, which contribute to the qqg final state. The corresponding insertion 
operator can be obtained by replacing uiq 
set {ttif}- 



in Eq. ( p.41| ), while retaining rriQ in the 



Explicitly, we obtain 



3(l,2,3|/3(e,/x2;fe,m,})|l,2,3)3 = x ^ 



2tt r(l -e) 



V2e^ 2e SQg J \ 2e^ 2e Sg^ / e eV6 3 / 
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QQ 
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1 rrin 
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1 2Cf — Ca 



^QQ 



l-Vr 



Np 
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'QQ 



rrir 



F=l Q 
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aux 
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where i/jk = Sjk/s and Qaux is an arbitrary scale that cancels against an equivalent term in 
G, as discussed after Eq. (|6.24|) . The function G contains the terms that do not diverge as 
either e — or rrip -C s, 



G{pi,p2,P3) = tCaIh 



4 



16. 



+ 3Cf In ■ 



2 2 



+ I (^Ca - Ir^Nf ] In 



+6Cf + -Ca - -T^Nf - Cftt^ - -Catt^ 
y y o 



6 3 

t^oo 6 



-Ca^Ih 
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Q 



In- 



1 



CA-ln 



In- 



'Qg 



Snn + 2?77,3i 



2Cf — Ca , 
H m 



\1 + 



In- 



+ ^Vf s)(sQ9, 0, niQ, {mp}; «:) + ^vf ^Hsq,, 0, mg, {m^}; «:). 
Combining Eqs. (pl37|) and ( p^ , we obtain 
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,(l,2,3|/3(e,/i2;fe,m,})|l,2,3)3+|A<3(Pl,P2,P3)| 



(l-loop) 
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27r 



:Ca 7£; - In 



47r/i^ 
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1 2Cf - C 
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2Cf — Ca 
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In 



- Ca In ^ Q - 7,; 



+i^finite(Pl,P2,P3)- (D.44) 

The integration of this expression over the phase space in Eq. (|D.8| ) with the observable 
Fj'^^ gives the three-parton cross section ci'^'^O't^l^ 

Comparing Eqs. ( p.44| ) and (p.40|) , we see that all logarithms of tuq cancel, leaving a 
cross section that is well-behaved for all masses. 
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If the observable to be calculated receives contributions from all final states, whether 
containing a heavy quark or not, it is also necessary to calculate the three-parton cross 
section for light quarks. Since our formalism guarantees a smooth small-mass limit, it is 
straightforward to set rriQ to zero in Eq. ( p.44| ), while keeping mp non-zero. 



D.3 ep eQQX 

Since the analytical formulae are more lengthy in this casep], we only give the formulae 
necessary to construct the auxiliary cross section and its integral. 

For simplicity, we start by assuming that the observable of interest requires the presence 
of a heavy quark pair in the final state. That is, we do not consider heavy-quark corrections 
to light-quark processes. Thus, at lowest order, there is only one process, e + g{pa) — > 
e -|- Q{pi) + Q{p2), with colour averaged matrix element 1/{N^ — l)\M.3{pi,p2, —Pa)\^- 
Note that this matrix element is not identical to the one given earlier for the process 
e~^e~ — > QQg, only because we averaged over the orientation of the three-jet system in 
the e^e~ rest frame. If the full angular information had been retained, then |A^3p could 
simply be obtained by crossing the e~^e~ annihilation result. The same comment applies 
to the tensor T^^, to be defined shortly. 

Three NLO real emission processes contribute: (a) e+g{pa) — > G + Q{Pi)+Q{P2)+g{P3)', 
(Jo) e + q{pa) ^ e + Q{pi) + Q{p2) + ^(ps); and (c) e + q{pa) ^ e + Q{pi) + Q{p2) + (/(pa)- 
The last is identical to the second so we do not explicitly consider it further. 

For process (a), we have to evaluate six dipole contributions, ^31,2, "^^32,1, "^^31, ^^32; T^f" 
and T)'^. The first two are constructed in exactly the same way as in the e+e~ processes 
already considered. The third dipole contribution, Dg^, with a final-state emitter and 
initial-state spectator, is given by 



with 



Pa-P3+Pa-Pl-P3-Pl m^«^ 

X31,a = ■ , (D.46) 

Pa-P3+Pa-Pl 

P'a = ^31,aP'a, P^l = P'3 + P'l ' ' ^^M. (D.47) 



and (y^^Q^) in Eq- ( ^.501 ). The dipole contribution "Dgg can be obtained from by the 
replacement pi ^ P2- 

The fifth dipole contribution for process (a), T>'f' ^ with an initial-state emitter and 
final-state spectator, is given by 

= ^ \ ^ X ^^T,^{P3UP2. -Pal (D.48) 

2pa-p3 2 Xzi,a 1 



ttA NLO Monte Carlo program, based on the subtraction method, to compute infrared and coUinear 
safe observables in this process is presented in Ref. pG|. 
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with x^i^a, Pa and p^i given by Eqs. (|D.46| ) and (|D.47| ) and {fj,\'Vf'^^\h') given in Eq. (|5.85|) . 
The tensor T^^ is the squared amphtude for the LO process not summed over gluon polar- 
ization, normahzed so that —g^^T^,^ 
from 'D'f' by the replacement pi ^ P2- 



\A4s\ . The dipole contribution can be obtained 



In process (b), we only have to evaluate two dipole contributions, T>1^ and 'Df'^ (recall 
that we assume for now that our observable requires heavy quarks in the final state. The 
other two dipoles for this process, T'12,3 and 'Df2, in which the heavy quark- antiquark pair 
is replaced by a massless gluon, do not therefore contribute). They are given by 



a3(b) 



1 



2pa-p3 2 X3i,o 



(/iivrii^) X 



-'T'^u{P3l,P2, -Pa), 



(D.49) 



with X3i^a, Pa and psi again given by Eqs. ( p.46|) and ( p.47|) and (/i|Vf'"'^ given in 
Eq. ( |5.83 ). The dipole contribution can be obtained from Vf^ by the replacement 

Pi ^P2- 

Next we need to evaluate the insertion operator /2+a(e, /x^, {pj, mj},pa), which cancels 
all the singularities of the one-loop cross section, and the finite operators related to the 
factorization of initial-state singularities, (^5 A^f! {Pi}^ xpa) and K2"' (x; {pi, mj}, xpa)- 

The expression for /2+a(e) is identical to that for Eq. (|D.42| ), except that {mp} 

is replaced by the empty set, {}, and k is replaced by 2/3. Note that our uniform notation 
(see the discussion after Eq. (C.27) of Ref. [|lT|) means that all dot products, Sja, remain 
positive and that s should be replaced by Q"^ = —{Pa — Pi — P2Y- 

The operators are independent of the presence of massive quarks in the final state 
and are given by 

2,a' (1, 2; XPa\PT [X] IJ.p; Pl,P2, XPa) 1 1, 2; XPa)2,a' = 

1 



a, 



—P''3^x)\n .. , 

27r X^SiaS2a - 1 



\M3{pi,P2, -XPa) 



2,a'(l, 2; XPa\PT [X] Hp] Pl,P2, XPa)\l, 2; XPa)2,a' 

1 



a. 



X 



1P99(^)1^ 

27r xyil^s^ " 



\M3{Pl,P2, -XPa) 



(D.50) 



(D.51) 



for incoming quarks, a = q, and gluons, a = g, respectively. 

The operators K"^'"" do depend on the quark mass and are factorization-scheme, but 
not -scale, dependent. They are given by 

Y 2,a' (1, 2; XPalKl'" (x, Pi, mg, ^2, mg, XPa) 1 1, 2; XPa)2,a' = 



27r 



K''(x)-K^i(x) + hQ 



2^F 

1 - x)si, 



c c 

-^ICl'^ix; Sia, mq) + 7^/C^'^(x; S2a, "^Q) 



+ iP/4(x) In ^^7" , + In 



'1 - X)s2a 



(1 - X)s2a + ml 
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x ^/_ ^ |.M3(Pi,P2,-m)P, (D-52) 

2,a' (1, 2; XPal^^f" (a;, PumQ,p2, rnq, XPa) 1 1, 2; XPa)2,a' = 



27r 



- 2mQ J sia + m?Q + 2m|, _ S2a - ^mq J S2a + m?Q + ^ml^ 




S2a 



m - 1 



(D.53) 



for incoming quarks and gluons respectively. The functions K"^"^ , /C^'^, /Cg'^ and P^gg are 
defined in Eqs. ( |6.56|) , (|6.59|) , ( |6.60|) and (|5.89| ) respectively, while K'^i^ define the factor- 
ization scheme [x) = in the MS scheme) as discussed in Ref. 



The ingredients given above are sufficient to construct a complete NLO calculation for 
any observable that requires the presence of heavy quarks in the final state. The dipole 
contributions are subtracted from the real matrix elements as in Eq. ( |6.6| ) to give a finite 3- 
parton integral. The operator /2+a(e) is inserted into Eq. ( |6.8D , cancelling the singularities 
in the one-loop matrix element to give a finite 2-parton integral. Finally, the operators 
(^j /^f) ^'^d K'^"' [x] are inserted into Eq. (|6.9| ) to give the finite remainder from 
initial-state factorization. 

However, for a complete description of all final states in deep inelastic scattering, one 
must also consider heavy-quark corrections to light-quark processes. Since we do not con- 
sider incoming massive partons, the only such correction comes from the 'QCD Compton 
process', e + g-^e + g + (7, followed by the gluon decay g ^ Q + Q. Rather than giving 
all the formulae relevant for a complete NLO calculation of the QCD Compton process, we 
assume that such a calculation already exists with Nf massless quark fiavours and give the 
additional terms that must be added to it owing to the presence of the Np massive quark 
fiavours. 

Firstly, we have to evaluate the two dipole corrections we neglected in process (b) earlier, 
1^12,3 and 'DI21 in which the massive quark-antiquark pair is replaced by a massless gluon. 
The first is constructed in exactly the same way as in the e+e" processes discussed earlier. 
The second, with a final-state emitter and initial-state spectator, is given by (recall that 
the momenta are labelled e + q{j)a) — ^ e + Q{pi) + Q{j)2) + ^(pa)), 

^12 = ^—TT^ \ ^ (/^iVg.Q.k) X ^T^.(P3, -Pa,Pl2). (D.54) 

2pi-p2 + 2mn 2 Xi2,a N, 



"Q 



c 



with 



_ Pa-Pl+Pa-P2-Pl-P2-ml 

Xl2,a — ; 5 (^JJ.OOj 

Pa- Pi +Pa-P2 

= Xl2,aP''a, = + " ( 1 " Xu,a)p''a, (D-56) 
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and (/i|VQ^Q_^ |z/) in Eq. ( 5.51 ). The matrix element for the Born process, e + q{pa) — »■ 
e + q{pi) + g{p2), not summed over gluon spin, is given by -^T^y{pi, -pa,P2)- 

Next we have to evaluate the additional contribution to the operator J(e) from the Np 
flavours of massive quark. 



2,a(l,2;a|5/2+a(e)|l,2;a)2,a 
= 2,a(l, 2; a\l2+a{e, {mF})\l, 2; a)2,a - 2,a(l, 2; a\l2+a{e, {})|1, 2; 0)2,0 

1 + PF,21 PF,21 /Q , 2 ^ 1 1 "^F 

^(3 + Pf,2i) - 9 111 — 

o ^ S21 



{Ar2i 
E 
F=l 



27r 



In- 
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E 
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1 + PF,2a 



2 

PF,2a 



. (3 + p|,2a)-^ln — 

3 ' 2 S2a 



X —\M3{pi,-Pa,P2) 
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(D.58) 



with ppjk = y^l — Amp/sjk and A'"^'^ defined to be the number of flavours for which sjk > 
Am^p. Note that this contribution is finite, so we have set e — > 0. It does however diverge as 
nip which is cancelled by a corresponding divergence from the virtual matrix element, 
yielding a two-parton integral that has a smooth small-mass limit. 

The operator P{x) is unaffected by the quark mass. Finally, therefore, we just have 
to evaluate the extra contribution to the operator K{x) from the Np flavours of massive 
quark, 

2,a'(l, 2; XpJ5K"''^'(x)|l, 2; XPa)2,a' 

a' 

= Yl 2,a' (1, 2; XPalK^'''' (X, {nip}) 1 1, 2; XPa)2,a' - E 2,a' (1, 2; XPalK'^'^' (a;, {}) 1 1, 2; XPa)2,a' 

(D.59) 



\f2a 



^x-TRj:^(<5(x+-x)-5(l-a;)) 



In 



1 + PF,2a 



1 / 1 - X + 2m?p/s2a 



S2a 



'1- 



^-(3 + PF,2a) - 9 111 — 



1 



1 — X 



-x)-pp2a r ^ Trl'^3(Pl,-a;Pa,P2)| 
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for a = q, with x+ = 1 — 4m|./s2a = p|'2a- Fo^' = 9i SK{x) is zero. Note that, as discussed 



at the ends of Sects. uM App. |B|, S2a is given by 2p2-pa, i-e. it is calculated 

from Pa, the momentum of the initial parton, rather than xpa, the momentum of the parton 
entering the Born cross section and from p2, the final-state momentum belonging to the 
boosted phase space. 

In the small-mass limit, since x+ 1, the first line of Eq. ( p. 60] ) does not contribute, 
so the whole contribution remains finite. It is straightforward to check that in this limit, 
it gives the same results as if the Np extra flavours were massless. 



These ingredients, together with the massless results in Ref. ||Tl[, are sufficient to provide 
a complete NLO calculation of all final states in deep inelastic scattering involving either 
massive or massless partons. However, as discussed in Sect. p.2| , care must be taken in taking 
the small-mass limit, since we do not include incoming massive partons. Logarithms of 
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will therefore remain in the final cross section, preventing a smooth small- mass limit. Since 
we have explicitly demonstrated that all the insertion operators do have smooth limits, 
this behaviour is isolated in the three-parton integral, J{da^ — da^). To achieve a smooth 
small-mass limit, the logarithmic behaviour of the three-parton integral can be evaluated 
and properly matched with a suitable definition of the heavy-quark parton distribution (see 
Ref. H). 
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